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Introduction 


In Unit 4 you met Newton's second law and the law of terrestrial gravitation. If we 
consider an object falling near the earth’s surface, and ignore forces other than 
gravity acting on the object, then these laws lead to a model of its motion in 
which the object’s acceleration is constant. If the height of such an object above 
ground level is x at time t then the mathematical statement of this model is 


dv _ 
ae 
ae ; ise ; 
where v = 5 is the velocity and g is the gravitational acceleration. This equation 
can also be written 
dx 
pie ene 1 


2 


rx; eae Le ee ; 
where ae is an abbreviation for the second derivative (eI of x with respect to 


t. Because Equation (1) contains a second-order derivative it is called a second- 
order differential equation, and as such is an example of the sort of equations that 
will be discussed in this unit. 


In Unit 4 we saw that Equation (1) can be solved rather simply, by direct 
integration. We have 


d |dx 
dt \dt 
Integrating with respect to t gives 


dx 


=-—-gt+A 


——— “<a 


where A is a constant. Integrating a second time gives 
x= —igt?7+At+B (2) 


where B is another constant. Equation (2) is the general solution of the second- 
order differential Equation (1). The constants A and B in the general solution may 
take any value—they are arbitrary constants. 


This example illustrates some general features of second-order differential 
equations. The appearance of two arbitrary constants in the general solution is 
typical of second-order differential equations. In Unit 2 you saw that one condition 
was needed to pick out a particular solution of a first-order differential equation. 
For a second-order equation, two conditions are required to pick out a particular 
solution. These conditions set values on the two arbitrary constants. In the 
example, we can find values for the constants A and B if we know the position x 


d 
and velocity = at time t = 0. Equation (2) can then be used to predict the 
position of the object at all future times t. 


The example also illustrates very well why second-order differential equations are 
important. They arise commonly in modelling. In particular, they arise frequently 
in dynamics, because dynamical models often involve the use of Newton’s second 
law. Mathematical statements which are derived using this law involve 
acceleration, the second time-derivative of position. However, our concern in this 
unit will be with the mathematics of solving second-order differential equations. 
Elsewhere in the course you will see many examples of their use in modelling. 


We have seen that the general solution of Equation (1) can be found by direct 
integration. However, this method only works for equations of the form 


d2 
ae PO 
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(where the integrals of f(t) and [Reva can be found explicitly). Most of this unit is 


devoted to describing a method of finding analytical solutions of a wider class of By ‘an analytical solution’ we 
second-order differential equations. The technique does not depend on integration mean a solution which can be 
(it is easier!). Rather it is a matter of recognizing that the equation is of a certain oe ee 
type, for which a set routine then provides a solution. 


The second-order differential equations with which we are mainly concerned are 
linear constant-coefficient equations. (I will define these terms shortly.) An example 
is 
d*y dy 
2— + 3— 
oe oe 
Equations of this type have the advantage of being relatively straightforward to 
solve analytically. 


—y =4sin 21. (3) 


Equations which are not of this special type may not be so cooperative. In such 
cases numerical methods may have to be used. The basis of one such method is 
described in Section 5. 


Terminology and notation 


You will need to become familiar with certain jargon associated with second-order 
differential equations. (You may recognize some of it from Units 1 and 2.) We 
shall be exclusively concerned with linear equations. 


A second-order differential equation is linear if it can be brought If we allow p to be the zero 
to the form 


function then Equation (4) 
includes the first-order linear 
d?y dy oe ney a 

x + g(x)— + r(x) y = f(x considered in Unit 2. In this 
P( V2 a ox ()y. =F) (4) unit we wish to confine our 
attention to second-order 
where p, g, r and f are known, continuous, functions, and p is not equations and so we impose 
the zero function. the condition that p is not the 
zero function. 


I shall make a habit of writing linear equations in the form of Equation (4), with Obviously, letters other than y 
all the terms involving the dependent variable (in this case y) on the left, so that ae x may also be used for 
the right-hand side is a function of the independent variable (in this case x) only. ae ee 


The equations we will be particularly concerned with are constant-coefficient 
equations, that is, equations for which p, q and r are constant functions. Equation 
(3) above is an example of a linear, constant coefficient, equation. 


We concern ourselves here only with linear constant-coefficient equations because, 
as I remarked earlier, these are relatively easy to solve. However, even when the 
coefficients are not constant the idea of a linear differential equation is still 
important because there are a number of general theorems about differential 
equations which apply to linear equations. 


There is one last piece of terminology which we shall need when we study methods 
of solution. An equation which can be brought to the form 


d’y dy 
p(x) a + q(x) = + r(x) y=0 


(that is, Equation (4) with the zero function on the right-hand side) is a 
homogeneous linear differential equation. If the right-hand side of Equation (4) is 
non-zero, then we have an inhomogeneous equation. 


Exercise 1 
(i) | Write down the order of each of the following differential equations. 
fy, —) dy 
—= — —+3=0 
(a) 7s (©) dx? re . 
dy 


dy dy 
—+3y=0 d) 2 = 3— 
(b) oe y (d) ae te me 
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(ii) | Which of the equations in (i) are linear? 


(ili) Which of the following second-order linear differential equations are constant- 
coefficient equations? 


d’y d’y 
(a) wat eye ae it 

d*y dy ; d’y dy 
b) 2—~ + 3— + 4y = sin d) 2— + cosx— + 3xy=e& 
©) dx? dx . (d) dx? - tay = 


(iv) Which of the equations in (iii) are homogeneous? 


[Solutions on p. 49 | 


Comment on theorems and proofs 


Certain important results about differential equations are stated in this unit as 

_ theorems. Since the objective of this unit is to teach you how to solve certain 
differential equations, proofs of theorems are not given. Theorems are stated where 
they give important information about the solutions of differential equations, and I 
will try to illustrate their importance. 


Many theorems require the condition that all the functions involved are twice 
differentiable throughout a suitable interval, and that this second derivative is 
continuous. I shall assume throughout the unit that all functions mentioned satisfy 
these conditions, without stating this explicitly in the theorems. 


Providing that you can solve the differential equations, understanding the 
theorems is only a secondary objective of this unit. No knowledge of the proofs is 
required. All the proofs omitted may be found in text books such as An 
Introduction to Linear Analysis by Kreider, Kuller, Ostberg and Perkins (Addison- 
Wesley 1966), and many also in M203: Introduction to Pure Mathematics. 


Study guide 


The first four sections of this unit are concerned with linear constant-coefficient 
equations. These are equations of the form 


d*y dy 
a—> + b— + cy = f(x 
TS aes I(x) , 
where a, b, and c are constants (and a ¥ 0). In Section 1 we examine the 
homogeneous case, where f is the zero function. The inhomogeneous case, where fis 
some non-zero function, is considered in Section 2. 


Section 3 is concerned with finding particular solutions that fit given conditions 
dy 
a 
this section. Section 4 is the television section, and examines the graphs of some of 
the functions that arise as solutions of linear constant-coefficient equations. Section 
5 concerns a numerical method for second-order differential equations. Section 6 
contains additional exercises. 


(such as y = 2 and 1 when x = 0). There is an audio-tape associated with 


The television section makes some reference to the material in Sections 1 to 3 and 
so it would help to have studied these sections before watching the programme. 
However, if you have not finished the first three sections by the time the television 
programme is due to be broadcast, make sure you have read the preliminary notes 
at the beginning of Section 4 before watching the programme. The audio-tape 

_ cannot be studied until you have completed Sections 1 and 2. Although there are 
five sections of expository text you should not find the unit excessively long since 
Sections 3, 4 and 5 are all fairly short. 
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1 Homogeneous equations 


1.0 Introduction 


In this section we are concerned with the solution of equations of the form 


The general 
d’y dy homogeneous 
a5 +b-+cy=0 (1) constant-coefficient 
dx dx second-order 
where a, b and c are constants, and a is not zero. I shall first describe (in differential equation. 


Subsection 1.1) the method whereby the general solution of an equation of this 
form can be found. Afterwards, in Subsection 1.2, I shall explain why the method 
works. 


1.1 The method of solution 


The auxiliary equation 
To indicate the idea of the method we consider the linear first-order differential 
equation 


One solution of this equation is y = e**. This suggests that we might try the same 
exponential form for the solution of the second-order Equation (1). 


To try y = e™ in Equation (1) we first calculate 


dy “= 
om = he 
and 
d’y 2 Ax 
ae 
so that 
iy 2a 
eee = cy = ai*e** + bhe* + ce” 


dx? dx 
| -= (ad? + DA + che. 


72038, 7 = "As is a solution of Equation (1) provided that 
al? + bA+c=0. 


This equation plays such an important role in solving the differential equation that 
we give it a special name: 


Definition 1 
The quadratic equation 


aj? +bA+c=0 


is called the auxiliary equation of the Differential equation (1). 


Example 1 
The auxiliary equation of the differential equation 
oy dy | 
ee ee Ee, 
ee dx Se 


347 -24+4=0. 
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Solving the differential equation 

You should recall that the solution of a quadratic equation with real coefficients 
(such as the auxiliary equation) may be divided into three cases: 

(i) |The equation has two distinct, real, solutions. 

(11) The equation has one real solution (two equal roots). 

(iii) The equation has two conjugate complex solutions. 

These three cases for the solution of the auxiliary equation correspond to three 
cases for the solution of the Differential Equation (1). In the procedure below I 
describe how to write down the solution of Equation (1) in each case. Later in this 


section (Subsection 1.2) I shall explain why the solutions take the form they do. 
For now, though, I am only concerned with the technique. 


Procedure 1.1 
To solve the homogeneous differential equation 


wf dy 
a aoa + ba + cy = 0 


where a, b and c are constants, and a ¥ 0. 
1. Write down the auxiliary equation 

ad? + bA+c=0 
and solve it. 


2. The solution of the differential equation takes a different form in 
each of three cases, depending on the type of solution of the 
auxiliary equation. 


(i) If the auxiliary equation has two, distinct, real solutions 
A =A, and 4 = d,, then the general solution of the 
differential equation is 


y = Ae** + Be*2*, 


If the auxiliary equation has equal roots, so that it has just 
one solution 1 = a, then the general solution of the 
differential equation is 


= (A + Bx)e™ 


If the auxiliary equation has complex roots, 4 = a + if and 
A = « — if, then the general solution of the differential 
equation is 


y = e*(A cos Bx + B sin Bx). 


In each case A and B are arbitrary constants. 


The remainder of this subsection consists of worked examples illustrating the 
procedure above. You may wish to read through these, or you may prefer to try 
the examples yourself first, and refer to their solutions only if you get stuck. 


Example 2 
Solve the differential equation 
2 
- = a + 2y = 0. 
Method 
The auxiliary equation is 
A* — 3A 4+2=0. 


Its solutions are 4 = 1, A = 2. This is Case (i) of Procedure 1.1 and so the general 
solution of the differential equation is 


y = Ae + Be™ 


where A and B are arbitrary constants. 
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Example 3 


Solve the differential equation 
2 


d*y 


oe ee 0. 


Method 
The auxiliary equation is 


A*#+4=0. 


The solutions of this are A = 2i, 4 = —2i. This is Case (iii) of the Procedure 1.1 
with « = 0, B = 2. Hence the general solution of the differential equation is 


3 


y =A cos 2x + B sin 2x You may wish to check the 


where A and B are arbitrary constants. solutions obtained in these 
examples by differentiation 
and substitution. 


Example 4 
Solve the differential equation 
d7u du 
4-— + 4— == (), 
Tee ae 


Method 
The letters used for the variables are different here, but this makes no difference to 
the method. The auxiliary equation is 


447+44+1=0. 


The a0 hand side is the perfect square (2A + 1)” and so the only solution is 
4 = —4. This is Case (ii) of Procedure 1.1, so the general solution of the 
diiconsial equation is 


u=(A + Bt)e~**, 


where A and B are arbitrary constants. 


Example 5 
Solve the differential equation 
ie a 


Method 
The auxiliary equation is 


247 — 31 +2=0. 
Using the formula for solving the quadratic equation gives 
34 9-16 3.2/7 
a 


Sak 


This is Case (iii) of Procedure 1.1, with « = -, B = = So the general solution of 


the differential equation is 


sens (44) « nn (7) 


ees 


¢=> esp 6) 


where A and B are arbitrary constants. 


Example 6 (Where 4 = 0 is a solution of the auxiliary equation.) 
Solve 
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Method 
The solutions of the auxiliary equation 4? + 3A = 0 are 1 = 0, A = —3. This is 
Case (i) of Procedure 1.1. The only difficulty lies in realizing that e°* = e° = 1, so 
that the general solution is 

y=A+Be*, 


where A and B are arbitrary constants. 


Exercise 1 (Practice in solving homogeneous differential equations.) 
Solve each of the differential equations below. 


d*y dy dy 
“ami cee eee = — + 2— 0 
(i) Ae “eo + 5y=0 (i v5 ai te +y= 

d*6 du 
Th eeeegres” e ea a ee 
(ii) 72 +90=0 (v) ie s cp + 8u= 

22 d 

(iii) wa 4z=0 (vi) 343% 0. 


[Solution on p. 49 | 


Sometimes we need to solve equations which have coefficients expressed in terms 
of unspecified constants rather than numbers. This does not affect the method of 
solution. 


Example 7 

Find the general solution of the differential equation 
d?x dx : 
esi ee =f} 
I + 2 Ox 


where q is a constant. 
Method 
The auxiliary equation is 
A? + 4m + w? = 0. 
The values of A satisfying this are 
| — 4a + ./16m? — 4m? 
4 == NV = (-2+ /3)o. 


We have distinct, real, roots (Case (i) of Procedure 1.1). The general solution of 
the differential equation is therefore 


wae =2+43 sat Saad o 
x = Ae 2+V3)0t 4 Bel-2-V/3)or 


where A and B are arbitrary constants. 


Exercise 2 

Find the general solution of each of the differential equations below. 
d? dy 

(i) aa t 20 + wy = 0 


d*x dx 
(il) SF z+ 60— + 50*x =0 
where is a constant. 


[Solution on p. 49 | 


1.2 Why the method works 


_We now look at why the expressions given in Procedure 1.1 for the general 
solution of a homogeneous equation (Equation (1)) are correct. Before looking at 
the specific expressions given for the solutions, I will discuss some ideas relevant to 
_ the solution of linear differential equations in general. 
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Definition 2 


If u, and u, are functions, then a linear combination of u, and u, 
is any function u of the form 


u(x) = c,u,(x) + C,u,(x) 


where c, and c, are constants. 


Example 8 
Suppose u,(x) = x?, u,(x) = sin x. The functions 
u(x) = 3x? + 2 sin x, 
u(x) = 2x? — 4 sin x, 
are linear combinations of u, and u,. On the other hand the functions 
u(x) = x? sin x, 
u(x) = x? + sin 2x, 
u(x) = 2sinx + x*, 


are not linear combinations of u, and u,. 


Theorem 1 


Suppose that the functions u, and u, are solutions of a This theorem is fairly easy to 
homogeneous linear differential equation. Then any linear sarge 

combination of u, and u, is also a solution of this differential 

equation. 


The next example illustrates how this theorem can help us validate the method in 
Procedure 1.1. 


Example 9 
Consider the differential equation 

u"(x) — 3u'(x) + 2u(x) = 0. (2) eel that u’(x) means 
Each of the functions u,(x) = e* and u,(x) = e”* is a solution of this differential A 


equation. (This follows from the fact that 4 = 1 and A = 2 are solutions of the 
auxiliary equation 4”? — 3A + 2 = 0.) We now know, from Theorem 1, that any 
linear combination of u, and wu, is also a solution of Equation (2). That is, any 
function 

u(x) = Au,(x) + Bu,(x) = Ae + Be?*, 


where A and B are constants, also satisfies Equation (2). For example, one 
solution is 2e* + 3e?*. 


In each of the three cases described in Procedure 1.1, the formula given for the 
solution of Equation (1) is a linear combination of two functions. 

Case (i): Ae** + Be*2*. 

Case (11): Ae** + Bxe™. 

Case (111): Ae*cos Bx + Be** sin Bx. 
Theorem 1 shows that so long as we know that the individual functions (e*'*, e*2* 
in Case (i); e and xe in Case (ii); e“sin Bx and e*“cos Bx in Case (iii)) are 
solutions of Equation (1), then these linear combinations must also be solutions. 


To check that the individual functions are in fact solutions is not too complicated. 
Let us look at this now. 


Checking the solutions 


The general homogeneous constant-coefficient linear second-order differential 
equation (1) may alternatively be written as 


au” (x) + bu’(x) + cu(x) = 0 (3) 
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(a notation which emphasizes that the solutions u(x) are functions). Now we saw 
in Subsection 1.1 that u(x) = e** is a solution of Equation (3) so long as A satisfies 
the auxiliary equation 

al? + bh +c=0. | (4) 
Let us now check that the functions specified in Procedure 1.1 are solutions of 
Equation (3) in the various cases. 
Case (i) 
In this case A = J, and A = J, are real distinct roots of Equation (4). Thus the 
functions u,(x) = e*'* and u,(x) = e*?* are already known to be solutions of 
Equation (3). 
Case (ii) 
In this case A = a is the only solution of Equation (4). Thus u,(x) = e* is known 
to satisfy Equation (3). The expression given in Procedure 1.1 for the general 
solution of Equation (3) also contains the function u,(x) = xe**. We must check 
that this is also a solution. We have 

u,(x) = xe™ 

u(x) = e* + axe™* = (1 + ax)e™* 

U3 (x) = ae™ + (1 + ax)ae™ = (20 + a?x)e™. 
Hence | 

au3(x) + bu3(x) + cu,(x) = e*(a(2a + a?x) + b(1 + ax) + cx) 

= e**((2ax + b) + (ax? + ba + c)x). (5) 

Now « is a solution of Equation (4); so aa? + ba +c = 0. 
Also, we know that the auxiliary equation has coincident roots, so that 


b 
b? — 4ac = 0. Hence its solution is « = — > and so 2aa +b =0. 


Thus the right-hand side of Equation (5) above is indeed ‘zero, and so u, is a 
solution of Equation (3). 


Case (iii) 
We know that A=a+i8 and A=a-—if are solutions of Equation (4) and 
so the complex functions 


u, (x) iad el + iB)x 
and 
uz(x) =e ™ 


are solutions of Equation (3). In this course we will normally restrict our attention 
to real solutions of differential equations, and so the complex solutions u, and u, 
are not suitable for our purposes. However, using Euler’s formula (see Unit 5, 
Section 4) we can write: 


1 
e*cos Bx = ae + —u,(x) 
z 2 
and. 
e*sin Bx = Le (x) — Lae (x) 
a bee 
Hence, by Theorem 1, the (real) functions e**cos Bx and e*sin Bx are also solutions 
of Equation (3). These are the functions appearing in the formula in Procedure 1.1 


The general solution 


So far we have seen why the functions given by Procedure 1.1 are solutions of 
Equation (3). To show that the procedure gives the general solution of Equation (3) 
we must show that there are no other (real) solutions. This is the consequence 

of another general theorem. 


To make this argument 
d 
precise we should define - 


when z is a complex variable. 
If w(t) = u(t) + iv(t), where 
u(t) and v(t) are real functions 
of t, then we define 


dw(t) _ du(t) | dv(t) 
a eee eS 


It is then possible to show 
that many of the usual rules 
of differentiation still apply. In 
particular 


de'Bt 
dt 


Furthermore, having defined 
the derivative of a complex 
variable, it becomes 
meaningful to talk about 
complex solutions of 
differential equations. 


= ie”. 


We shall normally invoke 
Theorem 1 in the context of 
real solutions, but the theorem 
is true for complex solutions 
as well. 
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The basic idea of the theorem is that once we have found two ‘genuinely different’ 
solutions of a linear second-order homogeneous differential equation, then the 
general solution is given by an arbitrary linear combination of them. To state the 
Theorem we need the following definition of what ‘genuinely different’ means. 


Definition 3 


Two functions u, and u, having the same domain are linearly 
independent if neither function is a constant multiple of the other. 


Example 10 


(i) If u,(x)=e* and u,(x)=e?* then u, and u, are linearly independent. 
For there is no constant A such that e?* = Ae* (or e* = Ae?*) for all real 
numbers x. 


(ii) If u,(x)=sinx and u,(x)=2sin x then u, and u, are not linearly 
independent since u(x) = 2u,(x). 


Theorem 2 


Suppose that D is an interval of real numbers and that the 
functions u, and u, both have domain D and are linearly 
independent solutions of the differential equation 


(for all x in D) (6) 


Suppose also that for every x in D, p(x) is not equal to zero. 


P(x)u"(x) + q(x)u'(x) + r(x)u(x) = 0 


Then the general solution of Equation (6) is Au, + Bu, where A 
and B are arbitrary constants. 


I have stated the theorem in a form applicable to general second-order 
homogeneous linear differential equations. Notice that a second-order equation 
with constant coefficients inevitably satisfies the condition p(x) # 0 for all x. 
(If p(x) = a = 0, Equation (6) is not a second-order equation!) 


My main aim in this section has been to explain why the method given in 
Subsection 1.1 for constant coefficient equations works. However, the theorems 
are informative for non-constant-coefficient equations also. 


Example 11 
The general solution of 
x*u"(x) + xu'(x) — u(x) =0 (x > 0) (7) 
is 
B 
u(x) = Ax + a (x > 0). 
Explanation 


Let u,(x) = x and u,(x) = ~ be functions which have the same domain (x > 0). 


We can show that these functions are solutions of Equation (7) by direct 
substitution. (I leave you to do this.) 


Now u, and u, are linearly independent functions. Also, Equation (7) is of the 
form described in Theorem 2 (for in this case D is the interval (x > 0) and so 
p(x) = x? is non-zero for all x in D. 


Hence, by Theorem 2, the general solution of Equation (7) is 
A 
u(x) = "4 + Bx (x >) 


where A and B are arbitrary constants. 


The function u, is a constant 
multiple of u, if both 
functions have the same 
domain and if there is a 
constant A such that 


U; (x) = Au,(x) 


for all x in their common 
domain. 


A definition of interval is given 
in the study comments for 
Unit 6 in Section 5 of the 
Preparatory Booklet. 
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Exercise 3 
Find the general solution of the differential equation 
dy dy 
2 —_ — xX— pao 
2x de oF Ee 0 (x > 0) 


Hint: first look for solutions of the form y=x" 


[Solution on p. 49 | 


Summary of Section 1 


The general solution of 


d*y _,dy 
aa + rr hey =0 
‘depends on the nature of the solution of its auxiliary equation 
al? +bA+c=0. 


The following table summarizes the method: 


General solution of differential 
equation. (In each case A and B 


Type of solution of auxiliary equation are arbitrary constants.) 
Two real roots: 4 =A, andd=A, y = Ae** + Be? 
Equal roots (one real solution): 2 = « y=(A+ Bx)e™ 


Complex roots: 4 = a + if and 1 = a — if y = e*(Acos Bx + Bsin Bx) 


This method may be justified by: 


1. verifying that the various functions (e*'*, e*?*; e, xe; e* cos Bx, e™ sin Bx) are 
indeed solutions in each case; 


2. referring to Theorems 1 and 2 stated in the text. 


The concepts of linear combination (Definition 2) and linear independence 
(Definition 3) for two functions were defined in the text. 


End of section exercises 


Exercise 4 
(i) | What is the auxiliary equation of the differential equation 
dy d’y 
3— —y-2—7, =0 
dx dx? 


(ii) Solve this auxiliary equation. 
(ii) Write down the general solution of the differential equation. 


[Solution on p. 49 | 


Exercise 5 

Find the general solution of each of the differential equations below. 
sa ay dy 
(i) Peg eee Best. ere ee SoS 

i Re ie 

(11) rm ae l6y =0 (iv) af 0. 


[Solution on p. 50] 


Exercise 6 
For which values of the constant k does the differential equation below have a general 
solution that involves sines and cosines? 
dy dy 
—; + 4k— + 4y =0. 
dx? eas 


[Solution on p. 50] 
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Exercise 7 


Given that y=x and y=vx? are solutions of 
2 


d*y dy 
hes = 
x 7 2x = + 2y=0 (for all real x) 


is it legitimate to use Theorem 2 to write down the general solution of this differential 
equation? 


[Solution on p. 50] 


2 Inhomogeneous equations 


2.0 Introduction 
In this section we study the solution of equations of the form 

au" (x) + bu’(x) + cu(x) = f(x) (1) 
where a, b and c are constants. These differ from the homogeneous equations 


studied in Section 1 in that a function f(x) has been introduced on the right-hand 
side. 


2.1 The method of solution: ‘particular solution 
plus complementary function’ 


In essence, the method of solution is similar to the method used in Unit 1 for 
solving inhomogeneous recurrence relations. The method depends on our being 
able to find a particular solution of Equation (1). Adding any such particular 
solution to the general solution of the homogeneous equation 


au"(x) + bu'(x) + cu(x) = 0 


(obtained by replacing f(x) in Equation (1) by zero) gives the general solution of 
Equation (1). 


Before examining the theoretical justification for this method I will work through 
an example that illustrates how it works in practice. 


Example 1 
Find the general solution of the equation 

u"(x) — 3u’(x) + 2u(x) = 6. (2) 
Method 


Step 1: We first solve the homogeneous equation obtained by replacing the right- 
hand side of Equation (2) by zero. That is, we find the general solution of 


u"(x) — 3u’(x) + 2u(x) = 0. 


Apart from a change of notation this equation is the same as the equation in 
Example 2 of Section 1. There we obtained the solution u = u,, where 


u,(x) = Ae* + Be*. 
(A and B are arbitrary constants.) 


Step 2: Next we find any one solution of the Inhomogeneous equation (2). This is 
not always easy, but in this case it is quite easily seen that the constant function 


u(x) = 3 


is a solution. (For if u,(x) = 3, then u(x) = 0 and u3(x) = 0. Thus 
u(x) — 3u,(x) + 2u,(x) = 6, as required.) 


Step 3: Finally we add the particular solution found in Step 2 to the general 
solution of the homogeneous equation found in Step 1. The result u,(x) + u,(x) is 
the general solution of the original Inhomogeneous equation (2). Thus 


u(x) = Ae* + Be* +3 


is the general solution of Equation (2). 
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The theoretical justification of this method relies on a theorem (Theorem 1 below). 
Although we are concerned with constant coefficient second-order equations in 
this unit, the method works much more generally, and so I will state the theorem 
for general second-order linear equations. To do so I will need the following 
definition. | 


Definition 1 
Let 


p(x) u(x) + q(x)u'(x) + r(x) u(x) = f(x) (3) 
be any inhomogeneous linear second-order differential equation. 


Then its associated homogeneous equation is 


p(x)u"(x) + q(x)u’(x) + r(x)u(x) = 0. 


Theorem 1 


Suppose that u, is any solution of Equation (3) and that u, is the 
general solution of the associated homogeneous equation. Then 
u, + u, is the general solution of Equation (3). 


We often need to refer to ‘the general solution of the associated homogeneous 
equation’ and so it is convenient to give it a special name: 


Definition 2 
The general solution of the associated homogeneous equation is 
known as the complementary function. 


With this definition we can make a succinct statement of the method for solving 
inhomogeneous equations of the form of Equation (1): 


Procedure 2.1 
To solve the inhomogeneous differential equation 


au"(x) + bu'(x) + cu(x) = f (x). 


-1. Find the complementary function using Procedure 1.1. 


2. Find a particular solution. 


3. General solution 
= Particular solution + Complementary function. 


Example 2 
Find the general solution of the equation 


u”(x) + 9u(x) = 9x. (4) 


Method 
The associated homogeneous equation is 


u"(x) + 9u(x) = 0 
which has the general solution u = u,, where 
u,(x) = A cos 3x + B sin 3x. 
This is the complementary function of Equation (4). 
A particular solution of Equation (4) is the function 
u,(x) = x. 


(For, if u,(x) = x, then u,(x) = 1, so u"(x) = 0. Thus u}(x) + 9u,(x) = 9x, as 
required. ) 
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The general solution of Equation (4) is, therefore, 
u(x) = x + Acos3x + Bsin 3x, 
where A and B are arbitrary constants. 


The finding of the complementary function is a question of applying the method 
from Section 1. However the finding of a particular solution is another matter. In 
the two examples above, I have carefully chosen equations for which a particular 
solution can be found ‘by inspection’. This sort of approach is not going to work 
in general, though. The next subsection looks at the finding of particular solutions 
in a more systematic way. 


Exercise 1 

Consider the differential equation 
d’y 
— + 4y =8. 
ae” 


(i) | Write down the associated homogeneous equation and the complementary function. 
(ii) Find a particular solution of the form y = constant. 
(iii) Write down the general solution of the differential equation. 


(Solution on p. 50] 


There is one point that sometimes causes confusion. An inhomogeneous equation 
has more than one particular solution (of course!). Does it matter which one we 

choose to add to the complementary function? It does not; the general solutions 

obtained may look different, but are actually equivalent. 


2.2 Finding a particular solution 


I shall not describe a general procedure for finding a particular solution of 
Equation (1). The form of the particular solution depends on the form of the 
function f(x) on the right-hand side of the equation. I shall describe the sort of 
functions that work for certain specific forms of f(x). On the whole, the approach 
is ‘try something similar to f(x). I will start with a fairly simple case. 


fislinear (f(x) =kx+ J) 

The first example I will consider has f(x) = 4x + 2. The basis of the method is to 
try a particular solution which is similar to f In this case fis linear so I will try an 
arbitrary linear function y = mx + n, and find (if possible) suitable values for the 
constants m and n for which this is a solution. 


Example 3 
Find a particular solution of the equation 
d*y dy 
a =4x+2 
3 s 2 Ps +y=4x+ 


Solution 
We try a solution of the form 


y=mx-+n 


where m and n are constants to be chosen so that the differential equation is 
satisfied. As a preliminary to substituting this function into the differential 


d d? 
equation we note that a m and eas = 0; hence 
dx dx 


2y od 
ee 0 nt On +0) 


ee 


= mx + (n — 2m). 


For the function mx + (n — 2m) to be equal to the function 4x + 2 on the right- 
hand side of the differential equation, we must have 


mx + (n — 2m) = 4x + 2 (for all x). 
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This can only be true if 

m=4 
and 

n— 2m = 2. 
We can solve these simultaneous equations (by substituting m = 4 into 
n — 2m = 2) to give m= 4, n= 10. Thus 

y= 4x + 10 

is a particular solution of the given differential equation. 


It is always sensible to check calculations such as these by substituting back into 
the differential equation as follows. 


Check: If y = 4x + 10 then 


d’y dy 
oat 2a FV HSK OW 2A 


= 4x + 2, 


as required. 


Exercise 2 


Find and check particular solutions of the form y = mx + n for each of the differential 
equations below 
dy 


es 
+ 2y = 2x +3. ws te 


d’y dy 
a = dx 


dx? “dx 
[Solution on p. 50] 


(i) 


(Notice that in part (ii) of the above exercise you still need to look for a solution 
of the form y = mx + n. Although the function f(x) is just a multiple of x in that 
case, there is no solution of the form y = mx.) 


f an exponential (/(x) = ke**) 
In the next example, f(x) = 3e?*. ‘Try something similar’ to f(x) again works. This 


time we try an exponential, y = me?*, and determine the value of the constant m. We try an exponential with 
the same exponent as f. 


Example 4 
Find a particular solution of 
d*y f+ 
Pa + y = 3e 
Solution 


We try a solution of the form 


_ y=me*, 
ay ee S 2 ee oe 5, 
Then 7 = 2me , and — a = 4me ae + me** = 5me**, 


Thus y = me** is a solution of the differential equation if Sme?* = 3e2*, that is if 


3 
m= 5 Hence 


is a particular solution of the given differential equation. 


Exercise 3 
Find a particular solution of the equation 
d*y dy 
2— — 2— = 2e~*, 
i eS 


[Solution on p. 50] 
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f sinusoidal (f(t) = kcosat + /sin@t) 


The most important case so far as applications are concerned arises when the 
function f is a sine or a cosine or a linear combination of them both. In the next 
example f is given by f(t) = 2sin3t. Although there is no cosine term we look for a 
particular solution of general sinusoidal type y = mcos 3t + nsin 3t. . 


Example 5 
Find a particular solution of 
d*y | dy 
—— + — + 2y = 2sin3t. | 
712 < + 2y eae (5) 
Solution 


This time we try 


y = mcos 3t + nsin 3t. 


Then 
dy = 3ncos 3t — 3msin 3t 
dt 
and 
d2 
nA = —9mcos 3t — 9nsin 3t 
dt 
SO 
d? d 
— cs ead + 2y = (—7m + 3n)cos 3t + (—3m — 7n)sin 3t. 
dt dt 
We require this to equal 2 sin 3t. So we need 
—7Tm + 3n=0 
— 3m — 7n = 2. 
Ee 3 7 eee 
These equations give m = Tae and n= ay So a particular solution is 
3 ae? 
y= Tag 3t — a 3t. 


Before asking you to do any examples, I will show you an alternative approach for 
this example. This uses complex numbers, and is easier so long as you are 
confident with their use. The method is based on the idea that any sinuosoidal 
function can be represented by a phasor (see Unit 5, Subsection 4.3). We look for 
a particular solution of Equation (5) which is a sinusoidal function (as in the 
solution above). This time, however, we do the calculations using phasors. 


Remember that the phasor of the sinusoidal function 
y = mcosat + nsin wt 

is the complex number z for which 
y = Re(ze’”) 

and is given by 


zZ=m-— in. See Unit 5, Subsection 4.3, 
| Example 4. 
With this in mind, we can look for a solution of Equation (5) of the form 


y = Re (ze*") (6) 
where z is a complex number to be chosen so that the differential equation is 
satisfied. We have 

dy 


a Re (3ize*"’) 
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and 
2 


d*y : 

dt = Re (— —9ze* g 3 

So (6) is a solution of Equation (5) so long as 
Re((—9 + 3i + 2)ze*") = 2sin 3t 


Now the phasor of 2sin 3t is —2i so if we write the right-hand side in phasor form 
we have 


Re((—9 + 3i + 2)ze*") = Re(—2ie*"). 
This holds if 
(—7 + 3i)z = —2i 
SO | 
4 
7-3 29 +29 


With this value of z, Equation (6) gives a particular solution of the Differential 
Equation (5). This particular solution is a sinusoidal function with frequency 3 and 


Z 


phasor as +h It can therefore be written 
7 
y=- — = 008 3t — 55 sin 3 (7) 


This alternative approach may look a little complicated at first but in practice we 
can abbreviate the calculation considerably. You will find that you can get from 
Equation (5) to Equation (7) fairly easily by applying the following general 
procedure: 


Procedure 2.2(a) 
To find a particular solution of 


pa dy 
ad + b—+cy=kcosat + Isinat. 


“de | dt 
1. Write the right-hand side in phasor form 
Re ((k — il)e'*). 


2. Try a solution in the form 


y = Re (ze). 
3. Substitute in the differential equation and simplify to obtain 
(—w7a + iwb + c)z =k — il, 
and solve for z. 


4. Then the particular solution is the sinusoid with frequency @ 
and phasor z. That is, if z = X + iY, the solution is 


y = X coswmx — Ysin wx. 


Although the exercises in this unit can be done using either of the methods I have 
described for Example 5, I recommend that you get used to the complex number 
method. It is always quicker, and sometimes much quicker (for example, for the 
sort of problems discussed in Section 4). 


Exercise 4 

Find a particular solution of the equation 
d’y dy 
ae + 3 = = 9cos 3t 


[Solution on p. 50 | 
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Exercise 5 (Optional, additional practice.) 
Find a particular solution of the equation 


d’y 
rs 4y = 8cos2t + 16sin 2t. 


[Solution on p. 51] 


In this subsection we have looked at the problem of finding a particular solution 
of Equation (1) for certain specific forms of f(x). We have seen that the following 
procedure usually works 


Procedure 2.2 

To find a particular solution of 
A + bY + cy = fix): 

1. If f(x) = kx + I, try 
y=mx +n. 

a Ni) = he™, try 
y = me™., 

3. If f(x) = kcoswx + Isin wx, try 


y = Mcosa@x + nsin wx, 


or replace the right-hand side by Re((k — il)e’®*) and try 
y = Re (ze’*”) 


as described in Procedure 2.2(a). 


I pointedly say that this procedure ‘usually’ works, for reasons we shall now 
examine. 


2.3 Exceptional cases for particular solutions 


Occasionally, the functions suggested in Procedure 2.2 fail to provide a particular 
solution. Let us see an example of this difficulty, and how it may be resolved. 


Example 6 

Find a particular solution of the equation 
dy 
oe ees 4 nee 2x 
a2 y =2e 

Solution 


According to Procedure 2.2 the obvious guess is y = me”. If we try this, 


d d? i 
then 2 = 2me* and 3 = 4me?*, giving 

ea aoe 8 

Pe — 4y = 4me — sme" = v. 


So there is no value of m which makes y = me”* a solution of the differential 
equation. What has gone wrong? The trouble is that the solution suggested in 


Procedure 2.2 contains the function e2* which is a solution of the associated 
2 


d 
homogeneous equation =~ — 4y = 0. It follows that me” is part of the 


complementary function and so cannot be a solution of the inhomogeneous 
equation. 
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When Procedure 2.2 fails in this way we can usually overcome the difficulty by 
multiplying the function suggested in the procedure by x. In this case we try 


y= me 
Then 

dy 2 2 2 

qe = me * + 2mxe** = m(1 + 2x)e?*, 
and 

d*y 2x 2x 2x 

—* 2me** + 2m(1 + 2x)e** = 4m(1 + x)e?*, 
so that now 

d*y : 2x 2x 2x 

— 4y = 4m(1 + x)e** — 4mxe** = 4me?*. 

es 
This can be made equal to 2e?* by choosing m = > which gives the particular 
solution 
1 
y = =xe”*, 


2 


In this example we see that the approach described in Procedure 2.2 does not 
work, but that by ‘multiplying by x’ we can still find a particular solution. 


Let us now look at two other examples where similar difficulties arise, and how 
they may be resolved. 


Example 7 
Find a particular solution of the equation 

d2 

7 + Oy = sin 3x. (8) 
Solution 


Following Procedure 2.2 we would try y = mcos 3x + nsin 3x. However, since 


both cos 3x and sin 3x are solutions of the associated homogeneous equation 
d*; : 

i + 9y = 0, this would just give zero when put into the left-hand side of 
Equation (8). To overcome this difficulty we again multiply the function suggested 


in Procedure 2.2 by x. That is, we try 
y = x(mcos 3x + nsin 3x). 

We can simplify the calculation by expressing this in the phasor form 
y = Re (xze**) 


where z = m — in is to be chosen so that the differential equation is satisfied. We 
have 


dy : 3ix 

= ee Re ((1 + 3ix)ze”*) 
and 

d*y | . 3ix 

Pe Re ((6i — 9x)ze”"*) 
SO 

d’y 


a2 + 9y = Re (6ize*”). 


The right-hand side of Equation (8) can be written Re(—ie**). If this is to be 
equal to Re(6ize*"*) we require 


6iz = —-i 
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that is 
ee 
=e 
Substituting this into y = Re(xze*”) gives the particular solution 
1 
= —- 3x. 
y gx cos 3x 
Example 8 
Find a particular solution of the equation 
d*y _ dy 
oe 5 aay m= 2x +2 (9) 
Solution 


This time we have a linear function on the right-hand side, so the first step is to 
try 


y=mx+n (10) 
dy d’y me 
as suggested in Procedure 2.2. If y = mx + n then —* m and oo os 0 giving 
d*y _ dy 
— + 2— = 2m. 
Be aga 


Thus, for (10) to be a solution of Equation (9), we require 2m = 2x + 2 (for all x). 


Unfortunately this is impossible since 2m is a constant where as 2x + 2 varies with 
x. This shows that the differential equation does not have any particular solutions 
of the form y = mx + n. (As with the two previous examples, Procedure 2.2 breaks 
down here because part of the function y = mx + nis a solution of the associated 


homogeneous equation. The part of the function which causes the trouble is y = n, 
2 


d 
for this is a solution of the associated homogeneous equation a3 + 2 = 0 and 


therefore does not survive substitution into the left-hand side of Equation (9).) To 
overcome this difficulty, I will again multiply the solution suggested in Procedure 
2.2 by x and try 


y = mx? + nx. 
d aoe 
Then = = 2mx + n and es 2m, so that now 
x 


dx? 
Soe a 
af + on = 2m + 2(2mx + n). 


For y = mx? + nx to be a solution we require 
2m + 2(2mx + n) = 2x + 2, 

that is 
4mx + (2m + 2n) = 2x + 2. 


This can only be true for all values of x if m and n satisfy the simultaneous 
equations 


4m = 2 
2m + 2n = 2. 


Hence m = 3 and n = 5. Thus the differential equation has the particular solution 
= = (x + x) 
y ae ,) . 


In each of the examples in this subsection the method described in Procedure 2.2 
did not work because the solution we tried contained a term that was a solution of 
the associated homogeneous equation. In each case we were able to find a solution 
using the procedure on the next page. 
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Procedure 2.3 
If the approach to finding a particular solution described in 


Procedure 2.2 fails to produce one, then try x times the function 
specified there. 


Exercise 6 
Find a particular solution of each of the differential equations below. 
ee | 
(i) aa + 4y = 2c0s 21. 
d’y dy 
ii) —> —3— + 2y = 4e*. 
(ii) Ae ae + 2y = 4e 


[Solution on p. 51 | 


_ That is very nearly the whole story. The next exercise completes it. Whether or not 
you try the exercise yourself, do look at its solution. 


Exercise 7 (Challenge) 
Find a particular solution of the equation 


[Solution on p. 51] 


2.4 The superposition principle 


There is one final result that is relevant to finding particular solutions. I shall state 
the theorem using a compact notation for the left-hand side of a general second- 
order linear differential equation. I will write 


d* d 
L(y) = p(x) 5 + a(x) = + r(x)y 


Theorem 2 
Suppose that 


L(y) = f(x) 


and 


L(y) = g(x) 


are any two linear inhomogeneous differential equations, with the 
same associated homogeneous equation L(y) = 0. Suppose also 
that y = u(x) is a solution of the first equation and y = v(x) is a 
solution of the second equation. 


Then, if a and b are constants, y = au(x) + bu(x) isa solution of 
the differential equation 


L(y) = af(x) + bg(x). 


This result is known as the principle of superposition. It states that if the right-hand 
side of our differential equation is a linear combination of functions for which we 
know particular solutions, then the same linear combination of the particular 
solutions is what we want. 


Example 9 
Find a particular solution of the equation 


2 


y mee 
a © Sy = 4sin 3x + 13. 


MST204 6.2 25 


Solution 
A particular solution of the equation 
d*y 
aad + 9y = sin 3x 
1 . 
isy= — ra cos 3x (see Example 7). A particular solution of 
d*y 
7a 9y=1 


gosta | | 
Ss y= 9 (Check it!) Then Theorem 2 tells us that 


= gee aes : 
oe 9° 


is a solution of the given differential equation. 


When the right-hand side of an equation is the sum of two parts we can either 
solve for the two parts separately, as in Example 9, or together, as illustrated by 
the following example. 


Example 10 

To find a particular solution of the equation 
d? d 
oa +34 +4 = 4e* — 7sin 21, 


we would try 


y = le** + mcos2t + nsin 2t 
(or y = le* + Re (ze"). 


Exercise 8 

Find particular solutions of each of the differential equations below. 
d? d 

(i) fe ~ 2 + y = 4e* — 3e* (Use the result of Exercise 7). 


dx 
(ii) a + 100x = 9cos 4t + 100. 


[Solution on p. 51 | 


2.5 General solutions 


In Subsection 2.1 we discussed how the general solution of an inhomogeneous 
equation of the form given in Equation (1) on page 15 can be found by adding a 
particular solution to the complementary function. Now that we have developed a 
method for finding particular solutions, you should be able to find general 
solutions for suitable functions f (x). 


It is sensible when finding a general solution to find the complementary function 
first. This will provide advance warning of the sort of ‘exceptional case’ discussed 
in Subsection 2.3. 


Exercise 9 (Consolidation of all of Section 2) 
Find the general solution of each of the differential equations below. 


3 dy d 
w 4 sos - 40 = 2 (iii) aa z 3 + 2y = 2e — Se 
@ dy 
(ii) er. 2 -Y =e? oe ng & 2 2 4 208 + 20%y = Ssin 200 + 0? h 


(where w and h are constants) 


[Solution on p. 52 | 
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Summary of Section 2 


1. The general solution of an inhomogeneous linear differential equation is found 
by adding any one particular solution to the complementary function 
(Theorem 1). The complementary function is the general solution of the 
associated homogeneous equation. = 


2. To find a particular solution of an inhomogeneous constant coefficient 
equation, we try a particular expression for y, and see whether it works. The 
form to try for y depends on the form of the function on the right-hand side of 
the equation. The table below gives expressions that usually work. 


f(x) 


kx +1 
ke* 


mx +n 
me** 
MCOS Wx + nsin a 


kcoswx + Isinwx or 


Re (ze!”*) 


—— 


3. In some cases there is no particular solution of the form suggested in the table. 
This occurs when the function to be tried contains a term that is a solution of 
the associated homogeneous equation. 


If the expression given in the table for y does not work, then try x times the 
given expression. 


4. If the right-hand side of the differential equation is a linear combination of 
functions given in the table above, then a particular solution may be found by 
trying the same linear combination of the expressions given there for y 
(Theorem 2). 


End of section exercises 


Exercise 10 


What form of expression for y should you use, in order to find a particular solution of each 
of the differential equations below? 


e d? ese d? one 
(i) a = 4y = *, (iii) a sit 4y =< = 
2 ee d? . . d? . 
(ii) 2 — 4y = sin 3x, (iv) re + 4y = sin 2x + 3x. 


[Solution on p. 53] 


Exercise 11 

Find the general solution of each of the differential equations below. 
ee dy 

Oa 

ee 

(i1) ro oe 4y = sin 3x. 


[Solution on p. 53] 
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3 Initial and boundary conditions (Tape Section) 


In this section we will look at solutions of second-order differential equations 
satisfying additional conditions. In Unit 2 you saw that the imposition of an extra 
condition (such as y = 2 when x = 0) on the solutions of a first-order differential 
equation usually picks out one particular solution by sétting a value on the one 
arbitrary constant in the general solution. 


Now the general solution of second-order differential equation has two arbitrary 
constants. So to pick out a particular solution from the general solution we might 
expect to require two additional conditions. This is indeed the case. 


In the first part of the tape I will work through two examples which illustrate how 
we can find solutions satisfying given initial conditions. These examples will also 
provide you with an opportunity to review the techniques in Sections 1 and 2. In 
the second part of the tape I will discuss the uniqueness of solutions Satisfying 
given conditions. 


Start the tape when you are ready. 
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© Excample | 


Find the particular solution of 


Pig 
SS... FEE ee 
ax doc 
Which satishes y=O ana 3... mal when SOQ 
cl oc 
Solution 


au sciliary 


equation [zx ~ Zt A-1t) =O 


Se or A= 2 


genera \ 
Solution 


=| when, se = 5 


so 
i2He «+ 262° 


O=A+B () 
1=>aArt2B (2) 


Parbcvlar 
Solution 


a Example 2 


Find the particular solution of the inhomogeneous equation 


for which U 3 when t = O 
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ee Inhomogen eous equati ons 


Find the complementary function. 
Find any particular solution. 
Add these to get the general solution. 


Then apply the conditions to obtain the required 
particular solution. 


Note: the particular SelutionS obtained in steps 2 and 4 
are nef necessarily the Same. 


Avoaliacy- 
equation is 
A*+G4=0 

Se fe = St 


associ ated 
homogeneous 
equation 


th erv 
dt 
and 
d*y 2 
t2 
So 


ay —— t 
ie a = £2. 2:6. 


So we require L=%o, m= YQ and nz0 giving 


particulter 
Solu Ciore 


4= toe + Mt 


U=Acos 3+Bsin3trietge 


general 
SolvClior. 


30 MST204 6.3 Tape 


| 5 Applying the conditions 


= 
General solution Y= Acos 32t+Bsin Bt+ hoe +Yat 


Y=O when t=o 


O =A cCosO+Bsino+ Ke + %xO 


=A -X% 


differentiate 


So the particular selution satisfying the conditions js 


Meet & & . 1 <€ 
BS 7 a4g SI Se (7G ces St tee + 6t 


—_.—_—Y 


from 
Complementary 


Function 


original 
particular 
Solutior 
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© Checking IuLion. 


PROBLEM: solve d*y + Ay = e +t ; y= 0, = =} when t<=oO 
TS : 
7! 2% < 
SOLUTION : MN * oy6 St ai Cnt, BE ee +4t 


Checking _the expression, satisfies the equation. 


is Sin 3t and cos 3t are solutions of 


a Y= lee + qt is a particular solution 


a t a oS 
So the expression satisfies the differential equation 


Check the Solution satishes the conditions 


: 


dy a Lu a5 
<2. SS CO + = 
When = 4 g 7 1 a. 


Now work through Exercise 1 and check your solutions to make sure that you have 
understood the ideas so far. 


Exercise 1 
In each case find the solution of the differential equation satisfying the given conditions. 


(i) u"(t)+4u(t)=0; u(n/2)=0 and = u'(n/2) = 1, 


d? d d 
(ii) aa ~ 35 + dy = 4 y= 4 and = = —1 when x =0 


[Solution on p. 53] 


Now try Exercise 2, which is an introduction to the next part of the tape. Its solution 
is covered in the tape commentary. 


Exercise 2 
Can you find a solution of 


u”(x) + 4u(x) = 0 
such that u(0) = 0, u(xz/2) = 1? 


After you have tried Exercise 2, restart the tape. 
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es Answer to Exercise 2 


The general Solution is U(X) =ASin 2% +B cos 2x 


O=A sino + B CcosO =A sin T+ B cos 
so 


Conditions Same Valve of =< 


uwu(a)=kK and u(b) = 


Boundary or ak diferent 
Conditions valueS of 


w(a)=K and uw (b) = 


© Theorem |: [Escistence and uniqueness theorem for 
initial Condition problems 


let 
p(x) uw (x) +q4()u (x) + r(eute) = F(¥) 


be any linear second-order differential equation, for 
which p(x) does not take the value 2ero for any valve 


of x. Then there iS one and only_one Solution of Chis 
aifferential equation satisfying the pair of initial Conditions: 


u(a) =K and uta} = 
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Exercise 3 

Identify each of the following as either an initial condition problem, or a boundary 
condition problem. Solve each problem. 

(i) u”(x) + 4u(x) = 0; u(0) = 1, u’(0) = 

(ii) u(x) + 4u(x) = 0; u(0) = 0, u(n/2) = 0. 

(iii) u”(x) + 4u(x)=0; u(0)=0, uw'(0)= 

(iv) u"(x) + 4u(x)=0; u(—2)=1,  u(n/4) = 


[Solution on p. 53] 


These examples bring up some points of interest. Boundary condition problems 
may be ‘well-behaved’, in that they produce a unique particular solution, as in (iv) 
above. Indeed, this is often the case. However, a boundary condition problem may 
have no solution, as we saw in the tape, or it may have infinitely many solutions, 
as in (ii) above. __ 


For a homogeneous equation, with initial conditions of the form u(a) = 0, u’(a) = 0 
there is no need to do any detailed calculation, as the solution must be the zero 
function (c.f. example (iii) above. This short cut cannot be applied to boundary 
condition problems, as there is no guarantee that the zero function is the only 
solution (c.f. example (ii) above). Nor can it be applied to inhomogeneous 
equations, for then the zero function is not a solution at all. 


Exercise 4 
Solve each of the following problems. 


d*y dy 

oe Te eet eerie pa and 2 =) a Xe. 
d? d 

(ii) a y =O and > =0 at x=0. 
d? dx d 

(iii) Ga t wa + wx = 0; ‘x =aand 7 =Oatr=0. 
where w and a are given constants. (You can use Exercise 2(i) of Section 1 here, if 
you wish.) 


[Solution on p. 54] 


Exercise 5 (Challenge—optional) 


(1) | There is an analogue of Theorem 1 for first-order equations. What would you expect 


it to say? 
(ii) | Can you find any solution(s) of 
dy 
—— 0 
"de = 


satisfying y = 0 when x = 0? 
(iii) How does the answer to (ii) relate to your ‘existence and uniqueness’ theorem in (i)? 


[Solution on p. 54] 


Summary of Section 3 


In this section we examined the problem of finding solutions which, in addition to 
satisfying a second-order differential equation, also satisfy two extra conditions. 


Such conditions are called initial conditions if they both occur at the same value of 
the independent variable. Otherwise they are called boundary conditions. 


For initial conditions we have an existence and uniqueness theorem which is 
stated in Frame 9. 


For boundary conditions there is no guarantee that a solution will exist. Even if a 
solution does exist, it is not necessarily unique. 
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End of section exercise 


Exercise 6 
Find the solutions (if any) of each of the following problems. 
(i) u(t) + 4u'(t) + Su(t) = 0; u(0) = 0, u'(0) = 2. 
(ii) u(t) + 9u(t) = 0; u(0) = 0, u'(n/3) = 1. 

d*y 2dy dy 


soa ae a ace y=0, Ci jane es 


State in each case whether it is an initial or a boundary condition problem. 


[Solution on p. 54] 


4 Oscillations and graphs (Television Section) 


4.0 Introduction 


In this section we study the graphs of certain functions that arise as solutions of 
second-order differential equations. We will concentrate on cases of importance in 
dynamics, particularly those where the solutions contain sinusoidal terms, and so 
are oscillatory. Many of the ideas discussed in the programme are also covered in 
the first three sections of the unit. Those ideas introduced in the programme, and 
not covered elsewhere in the unit, are summarized in Subsection 4.2. 


4.1 Preliminary work 
The programme divides into three parts. 
1. We study the differential equation 


d2 

=z ~ ky = 0. i (1) 
If k is positive the solution of Equation (1) involves exponential functions, but if k 
is negative it is sinusoidal. The programme reviews some ideas associated with 
sinusoids that were introduced in Unit 5. It also examines the relation of the 
amplitude and phase of the solution to initial conditions specified on the solution 
and its derivative. 


2. We study the homogeneous equation 


d? d 
oa _ 2a + w*y =0. (2) 

(The algebraic forms of the coefficients are chosen because they provide simplicity 

later.) We know that the solution of Equation (2) is oscillatory (in fact purely 

_ sinusoidal) if « = 0, because we then have Equation (1) with k negative. For 

0.< « < @ the solution is an oscillation whose amplitude decreases with increasing 

x, but if « > @, the solution is not oscillatory. 


Either way, so long as a > 0, the solutions decrease to zero as x increases. 
3: We study an inhomogeneous equation: 
d*y ay 
ie = vk + Sy = sin 2x. (3) 


Any solution of Equation (3) is the sum of two parts. One part is sinusoidal with 
the same period as sin 2x, and is the same for all solutions of the equation. The 
other part is a solution of 


This part is an oscillation whose amplitude decreases as x increases, and becomes 
very small as x becomes large. (This second part of the solution may therefore be 
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called transient.) For large x only the first part of the solution remains (and so this 
part is sometimes called the steady-state solution). 


To gain maximum benefit from the programme you will need to have studied the 
preceding sections of the unit. If you have not, you will find many of the ideas 
from these sections reviewed in the programme, but you will probably find that 
the pace of the programme is too fast for you to follow the details. Either way, 
remember that the important things to concentrate on are the shapes of the graphs 
of the solutions of Equations (1) and (2), and the way a solution of Equation (3) 
can be divided into two parts. 


The programme makes some use of the ideas concerning sinusoidal oscillations 
and phasors that where introduced in Unit 5. Exercise 1 serves to remind you of 
these. Before viewing the programme, work through this, and, if you have studied 
Section 1 of the unit, Exercise 2 as well. 


Exercise 1 
Find the phasor of the sinusoidal function 


f(t) = 4cos St — 3 sin St, 
and hence find the amplitude and phase of this oscillation. 
[Solution on p. 55] 


Exercise 2 

Solve the auxiliary equation of the homogeneous differential equation 
d’y y 
Ba t 2g, + ory =0. 


For what range of values of « and w are the roots real? 


[Solution on p. 55] 


As soon as possible after viewing the programme, read the notes in Subsection 4.2 
and then work through Exercises 3, 4 and 5. The remaining exercises may be left 
till later if you wish. If you have not studied the earlier parts of the unit it is still 
advisable to study the notes and try Exercises 3, 4 and 5 after viewing the 
programme, although there will be points that you will need to return to after you 
have studied the first three sections of the unit. 


Now watch the television programme ‘Good vibrations’. 


4.2 Good vibrations 


This subsection summarizes the new material in the television programme that you 
will be expected to know, and also takes some of the ideas a little further. 


1. An equation that will be useful in Unit 7 is 
2 


y 
a2 + wy == 0. (4) 
The general solution of this equation (see Section 1) is 
y = Ccoswx + Dsinawx, 


where C and D are arbitrary constants. You saw in Unit 5, Subsection 4.3, that 
this can also be written as 
y = Acos (wx + @). 7 (5) 


Thus (5), with A and ¢ arbitrary constants, is an alternative way of writing the 
general solution of Equation (4). This solution is a sinusoidal oscillation (see 
Figure 1). The constant in (5) is called the angular frequency of the sinusoid, and 
is determined by the differential equation itself. The amplitude A and phase @ of 


the oscillation depend on the initial conditions (the values of y and ed at x = 0). 


TV6 
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: A 
; Starting point 
' (\ depends on phased 
and amplituge A 
- 
amplitude 
ae. 


period 22 


Figure 1. The graph of the sinusoidal oscillation y = A cos (wt + ?) 


Fitting different initial conditions will have the effect of expanding or contracting 
vertically the graph in Figure 1 (varying A), and shifting the curve sideways 
(varying #) as in Figure 2. 


Figure 2. The effect of varying the constants A or @ in Figure 1 

2. An equation that will be useful in Unit 8 is 
dy ay 
— + 2a + wy = 0. 6 
dx? dx : (6) 


The nature of the solutions of this equation depends on the signs of « and 

a? — w?. Of particular importance for dynamical applications is the case when «@ is 
positive. The fundamental nature of the solution in four possible cases is 
summarized below. 


positive negative 


See the solution to Exercise 12 
at the end of this section for an 
Increasing explanation of the entries in this 


exponential table. 


Decreasing 


Ositiv : 
P : exponential 


Oscillation 
with increasing 
amplitude 


Oscillation 
with decreasing 
amplitude 


negative 


The graphs of typical solutions in each of these four cases are shown in Figure 3. 
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(a) a > 0, a? — w*>0 (b) a <0, a? -—w?>0 


(c) a>0, a7 —w? <0 (d) «<0, a?7-—@w?<0 


Figure 3. Graphs of various solutions of Equation (6) 


Let us concentrate for a moment on the case when a > 0, «2 — w?2 >0 

(Figure 3(a)). In this case the general solution of Equation (6) is a linear 
combination of two decreasing exponentials. By choosing different linear 
combinations we can find solutions fitting different initial conditions. These 
solutions have graphs which are variants of the graph in Figure 3(a). Figure 4 
shows the sort of graphs that we can obtain in this case. All the solutions have the 
property of tending to zero as x becomes large. 


Ya y 
x 
Jy J 
a x =, 
Figure 4. Graphs of possible solutions of Equation (6) when « > 0 and «? — w? >0, fitting 
various initial conditions. 


In the cases where « < 0, a” — w” > 0, we get solutions showing similar variations 
on the basic increasing exponential shape shown in Figure 3(b). In the oscillatory 
cases (Figure 3(c) and 3(d)), the general solution of Equation (6) can be written in 
the form y = Ae” “cos (vx + @) where v = ./w? — «? and A and ¢ are arbitrary 
constants. Changing the values of A and @ will alter the initial conditions 

satisfied by the solution, but the solution curve will have the same basic shape. 


38 MST204 6.4 


Figure 5. In oscillatory cases, changing the initial conditions changes the values of the 
arbitrary constants A and ¢. The period remains the same. 


3. We saw in the programme how a solution of 


os d . 
3 + = + Sy = sin 2x : (7) 


is the sum of two parts (see Figure 6). 


A 


(a) A solution of Equation (7) 


y 


(b) Its transient part 


(c) Its steady state part 


Figure 6. The solution (a) is the sum of two parts, shown in (b) and (c) 
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One part is a decreasing oscillation (Figure 6(b)), the other is a pure sinusoid 
(Figure 6(c)). For sufficiently large x the solution becomes effectively equal to the 
sinusoidal part. Whatever initial conditions we apply to Equation (7), this 
sinusoidal part of the solution is exactly the same. (That is, its amplitude, phase 
and frequency do not depend on the initial conditions.) The initial conditions 
affect only the decreasing part of the solution. 


The decreasing part of the solution is called a transient. The sinusoid to which the 
solution settles when x is large is called the steady-state solution of Equation (7). 


These ideas apply more generally. The general solution of 


d? d 
sts + 2a + w*y = kcos vx + Isin vx (8) 
dx dx 

is the sum of the complementary function—the solution of 
ry. ay 
— + 2a— 4y = 0, 9 
el + 20 om + w*y (9) 


and any one particular solution of Equation (8) (see Section 2). We have seen that 
if « > 0 the solutions of Equation (9) tend to zero as x becomes large. In other 
words, if « > 0 the contribution from the complementary function is a transient. 

_ There is one sinusoidal particular solution of Equation (8), and this is the steady- 
state solution. Any other solution of Equation (8) is the sum of a transient (from 
the complementary function) and this steady-state solution. The steady-state 
solution does not depend on the initial conditions. Only the transient term is 
affected by the initial conditions. 


4. Since the steady-state solution of Equation (8) (with « > 0) is a sinusoidal 
function, it can be represented by a phasor. The easiest way to find the steady- 
state solution is to calculate this phasor directly. The procedure for doing this is 
discussed in Subsection 2.3, and also in Example 1 below. Example 1 also 
illustrates how the phasor can be used to find the amplitude and phase of the 
steady-state solution. The angular frequency of the steady-state solution of 
Equation (8) is v, the angular frequency of the sinusoid on the right-hand side of 
the equation. 


Example 1 
Find the amplitude, phase and angular frequency of the steady-state solution of 
d*y dy 25 
a + *7 + 100y = 7 kX 


Use the results to write down the steady-state solution. 


Solution 
The steady-state solution is a sinusoid with the same angular frequency as the 
sinusoid on the right-hand side of the equation, that is 


25 
angular frequency = rv. 


Pe ae 
We are therefore looking for a solution of the form C cos oe + Dsin 4° or 


25 ; | 
equivalently A cos (7% + @). Since we are interested in the amplitude and phase 


we want the solution in the second form here. To do this we shall find the phasor 
z of this sinusoid by trying y = Re(ze°°/**), which leads (by Procedure 2.2(a)) to 


HE) 9) 


Solving for z gives 


| 
z= —5.(1 +i). 
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The amplitude and phase of the oscillation are given by the amplitude and phase 
of this phasor, that is 
2/2 


2D 
—3n 
4 


The steady-state solution is therefore 


y = 2v 2005 cl , 


4 4 


amplitude = |z| = 


phase = Arg z = 


Exercises 3, 4 and 5 are to be worked directly after watching the programme. 


Exercise 3 

Describe (in one sentence) the graph of any solution of the differential equation 
d7y 
a3 + wy =0. 


[Solution on p. 55] 


Exercise 4 
(i) What conditions must a and o satisfy if the genetal solution of the equation 


d’y dy 
pe Sd 2a + w*y =0 
| is to be a decaying oscillation? 
(ii) Sketch the graphs of the solutions of each of the differential equations below, that 
satisfy the specified conditions. 
dy dy dy 
= t35-+9y= 0; y= 


Py dy dy 
=e 0; y= oan ee 


(Do not solve the equations. Your sketches should show (I) behaviour near t = 0, 
(II) behaviour for large t, (Hu ) whether or not the solution is oscillatory.) 


(a) 


(b) 


[Solution on p. 55] 


Exercise 5 

Describe the solutions of the equation 
d? dy 
<3 + a= + wy = cos vx, 


where 0 < « < . In particular, comment on the long-term behaviour of the solution, and 
the nature of the part of the solution that is dependent on the initial conditions. 


[Solution on p. 55 ) 


The following exercises are to be worked after viewing the pi aii and after 
working the first four sections of the unit. 


Exercise 6 
Find the amplitude, phase and angular frequency of the particular solution of 


d’y 


Pe ei as 0 


d 
satisfying y = 4 and + = 9 at x = 0. 
[Solution on p. 55] 
Exercise 7 
Consider the differential equation 


u"(x) + w7u(x) = 
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2 
with the boundary conditions u(a)=k and u (« + i = |. Do these two conditions 
guarantee a unique solution? (Consider the cases k = | and k # | separately.) 


[Solution on p. 55] 


Exercise 8 
Find the particular solution of 
d?x dx 
res + 4m— It + w?x =0 
Se os : dx 
satisfying the conditions x = a and = Oatt=0. 


(Use the general solution of this equation found in Example 7 of Section 1.) 
[Solution on p. 55] 
Exercise 9 


Find the particular solution of 
2y 


a + peo ~ + wx =0, 


d 
where 0 <B <1 and 0 <o, that satisfies x = a and — = 0 at t=0 


Sketch the graph of this solution. 
[Solution on p. 56 | 


Exercise 10 


Find the amplitude, phase and angular frequency of the steady-state solution of 


d*x dx 
ae +o + 100x = 9cos St. 


[Solution on p. 56] 

Exercise 11 

Find the phasor of the steady-state solution of 
d?x 
dt ery 


where B and m are positive constants, and f < 1. Find the amplitude of this solution. 


+ 2BoX + w?x = cos vt 


[Solution on p. 56) 


Exercise 12 (Optional) 
Explain why the graphs of the solutions of the differential equation 


et. 

a + a +w*y=0 
take the forms specified in the table on p. 36. 
[Solution on p. 56] 
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Compare this with Theorem 1 
of Section 3. 
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5 A numerical method 


5.0 futrodeciion 


All the second-order differential equations you have encountered in this unit so far 
have been amenable to analytical solution, in that we have been able to express 
the solution in terms of a formula. Unfortunately, not all second-order differential 
equations which occur in practice can be solved analytically. For example the 
differential equation 


which arises in the study of the large oscillations of a pendulum, has no solution 
in terms of the elementary functions such as trigonometric or exponential 
functions. 


In such a situation we may resort to using a numerical method. In Unit 2 we used! 
Euler’s method to solve first-order differential equations of the form 


dy 
ax sera m(x, y ). 


This involves using a recurrence relation to generate values y, that give 
approximations to y(x,), where the difference x,, , — x, between successive values 
of x is a constant h (the step length). In Euler’s method we approximate the 


He ”* to obtain the recurrence relation 


derivative of y at x = x, by 
Vr+1 5 yr + hm(x,, y,)- 


In this section we will show how an extension of this method can be used to solve 
second-order differential equations. 


5.1 Euler’s method 


There are several methods of obtaining numerical approximations to the solution 
of a second-order differential equation. The method I will describe here is to 
replace the second-order equation by two first-order differential equations, and to 
apply Euler’s method to each of these. (Later in the course, we will see other, and 
better, methods than Euler’s for first-order equations. Use of these, together with 
the basic idea of replacing a second-order equation by two first-order ones, will 
provide more satisfactory methods for second-order equations.) 


The key idea is to introduce another variable z to Tepresent the first derivative of 
y. That is, we write 


—- 
dx 
If we now differentiate z, we obtain 
dz _ d’y 
dx x 


To see how this idea is used we look at the following example. 


Example 1 
Replace the second-order differential equation 
| dy dy 
ia 
by a pair of first-order equations of the form 
dy 


dx 
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and 
dz ; ; ae 
a (some function of x, y and z not involving their derivatives) 
x 
Solution 
We can rewrite the differential equation as 
dz dy 
— — 5— + 6y = 0. 
dx dx 26 


d 
Also, replacing = by z, we obtain 


d 
he + by = 0. 
dx 


Then the pair of first-order equations (written in a form suitable for applying 
Euler’s method) is © 


dy = 
dx 
dz 
dx 


Similarly any linear differential equation of the form 


= 5z — 6y. 


d’y 
pox 5+ ax) + r(x)y = f(x) 


(for which p(x) is non-zero for all x) can be written in the form 


oe 
aS (1) 
= = m(x, y, Z) 


where m(x, y, z) is some function of the three variables x, y and z not involving their 
derivatives. Many non-linear equations can also be treated in this way. 


Exercise 1 


Replace each of the following second-order differential equations by a pair of first-order 
equations of the form given in Equation (1). 


ee 
Ot 


(ii) + 32siny =0 


dt ed 
[Solution on p. 57] 


Having changed our second-order differential equation into a pair of first-order 
equations, we can set up the equations for Euler’s Method. 


Example 2 
Formulate Euler’s method for solving the equation 
d*y dy 


Pie Te + 6y = 0. | (2) 


Solution 
From Example 1 we can replace the equation by the two first-order equations 


dy 
=z | (3) 


—~ == 5z — Gy. (4) 
dx 
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To apply Euler’s method to Equation (3), we approximate the derivative of y at 
x = x, by 
eed Sp 
soar eee. 


Here y, is Euler’s approximation to y(x,), the value of the true solution at x,. The 
values of x, are given by 


(3) 


X, =Xq + rh. 
d 
Replacing in Equation (3) by the Expression (5) gives the recurrence relation 


Vr+1 = J, + hz, 

where z, is the approximation to z(x,). 

Treating Equation (4) in a similar way we obtain 
Z.4, = 2, + h(5z, — 6y,). 


Thus Euler’s method eeu: to the Second-order equation (2) yields the pair of 
recurrence relations: 


y r+i — y r oe hz, 

Z,41 = —6hy, + (1 + S5h)z,. 

Exercise 2 

Derive the recurrence relations for Euler’s method applied to the differential equations 
dy _ 

(i) 


dx 
dz 
ax —4z — 5y + cosx 


(ii) —=2 
— = —32siny 
x 


[Solutions on p. 57 | 


We are now almost ready to carry out some computations. All that remains is to 
see how the initial conditions fit into this scheme. Although my main purpose in 
developing a numerical method is to be able to attack equations that cannot be 
solved analytically, I shall deliberately illustrate the method with an equation that 
we can solve analytically. We will then be able to see how accurate the method is. 


Example 3 


Use Euler’s Method to obtain the value of y when x = 1 using a step length 
h = 0.1 for the differential equation 


d’y dy 
wo 0 


; ae a d 
with the initial conditions y = 1 and = = + when x = 0; 


Method 
We saw in Example 2 that Euler’s method gives rise to the recurrence relations 


Yr+1 fey yy - hz, 

2.41 = —6hy, + (1 + Sh)z,. 
The initial conditions can be written as y, = 1 and z, = 1. Thus with h = 0.1 we 
have 


Vr+1 =; yy +: 0.12, 


«= | and z, = 1). 
Ob + US2 Wo ee 


Zr+1 
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We can now use the recurrence relations to calculate y, and z,, then y, and z,, 
and so on. 


The following table (in which the results have been calculated to 3 decimal places) 
shows the results of the calculation. 


0.6y, + 1.5z, 


5.2 Accuracy of the method 
For the equation 


d’y dy dy 
dx? -2 7 ee 0; y= 1 and — = 1 when x =0 


considered in the previous subsection we can compare our numerical results with 
the analytical solution, which is 


y= 2 e2* = 
The following graph gives both the analytical and numerical results. 
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Figure 1 


These results are not particularly good for this value of h though it can be seen 
that the numerical results do follow the general shape of the solution. With smaller 
step sizes we would expect the results to be much closer to the true solution 
(provided we work to enough decimal places). However, the smaller the step-size 
the more work we have to do to evaluate y at x = 1. Thus in general we would 
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wish to use a computer to obtain our results. The following table gives the 
- approximation to y(1) obtained for various step sizes h. 


The correct value for y at x = 1 is 2e? — e? = —5.307. 


Exercise 3 


(i) | Continue the table below to find an approximation to y(3.2), where y is the solution 
of the differential equation 

d’y 

oe 


d 
satisfying the initial conditions y = 0 and = = 1 when x = 0. 


| - a 2 . 


ee 1.0 0.2 1.0 
0.2 0.2 1.0 | 04 0.96 
0.4 0.4 0.96 0.592 0.88 
0.6 0.592 0.88 0.768 0.762 
0.8 0.768 0.762 | 0.920 0.608 
1.0 0.920, 0.608 | 1.042 0.424 
1.2 1.042 0.424 1.127 0.216 
1.4 1.127 0.216 1.170 — 0.009 
1.6 1.170 | —0.009 | 1.168 —0.243 
1.8 1.168 | —0243. | 1419 | -O477 
2.0 1119 | -040 - ie ~0.701 
2.2 1.024 | —0.701 | 0.884 — 0,906 
2.4 0.884 | —0.906 | 0.703 — 1.083 
2.6 0.703. | —1083 0.486 

2.8 0.486 | 

3.0 

3.2 


(ii) | What is the exact solution of this equation satisfying the given initial conditions? 


(iii) Compare the results of Euler’s method, given in the solution to part (i) above, with 
the true solution. 


[Solution on p. 57] 


Although the results in Exercise 3 are not very accurate, the values obtained for y, 
do follow the general pattern of the solution. It would be reassuring if this were 
always the case, but unfortunately this is not so. This can be seen in the next 
exercise. 


Exercise 4 

Apply Euler’s method to find y(1) using a step length h = 0.2 for the differential equation 
d*y _,dy 

—~ + 12— + 20y = 0 

dx? dx . 


d 
with y = 1 and = = 1 when x = 0. Compare your results with the true solution. 
[Solution on p. 57] 


The reason for the behaviour shown by the application of Euler’s method in 
Exercise 4 is discussed in the unit on numerical solutions of differential equations. 
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Summary of Section 5 


Euler’s method for obtaining numerical approximations to the solution of a 
second-order differential equation is to put 


dy 
zZ=— 
dx 
and replace the second-order differential equation by a pair of first-order 


equations. Euler’s method as described in Unit 2 can then be applied to each of 
these first-order equations. 


End of section exercise 


Exercise 5 


Use Euler’s method with step length h = 0.2 to find an approximation to y(1), where y 
satisfies 


d 
with y = 1 and . = 1 at x = 0. Compare your results with a graph of the analytic solution. 
x 


[Solution on p. 58] 


6 End of unit exercises 


Exercises 1, 2 and 3 provide further practice in the main skill that this unit aims to 
teach—the solution of second-order constant-coefficient differential equations. You 
should try to make sure that you are confident on this point. If you are, you may 
like to try some of the other exercises in this section, which use the material in the 
unit, but do not concentrate solely on this central skill. 


Exercise 1 
For each of the equations below find the solution satisfying the given initial conditions. 


dy d d 
(i) aa +6 + 5) =0; with y= 2 and © = 1 atx =0. 
yg d 
(ii) wat 2 ae, with y =O and = = 1 at x =0. 
2 d 
(iii) T} — oy =3sin3r + 18; with y = 0, and = = 0 at ¢=0. 
0 7 do 
(iv) a (u > 0); ae = Land | = 2atu= 4 


[Solution on p. 58] 


Exercise 2 

Find the general solution of 
2 

(i) 2 + 16y = 4cos 4t + 8sin 4t; 
x 


oe oe : 
(ii) a + 4a, + 4Y = 3 t+ 8t. 


[Solution on p. 59] 


Exercise 3 


For each of the following differential equations sketch the graph of the solution Satisfying 
the given initial conditions. 


a 2 d 
(i) sate ty = 0; with y= —1 and = =—2 at x=0 
i a d 
i ee ee i ae a 
dx dx dx 


[Solution on p. 59 | 
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Exercise 4 


Suppose that x denotes the position of an object at time t, and that x satisfies the 
differential equation 


a oe 
ao ee he oe 


(i) | Sketch a graph showing the variation of position with time, a long time after the 
object is set in motion. Does the value of x ever become negative during this part of 
the motion? 


(ii) | Can you make more precise the meaning of the phrase ‘a long time after the object is 
set in motion’ in part (i). 


[Solution on p. 59] 


Exercise 5 


For what values of the constant 4 can a solution other than the zero function be found for 
the boundary condition problem 


u”(x) + A?n7u(x)=0; u(0)=0, u(l)=0? 
[Solution on p. 60] 


Exercise 6 


This exercise looks at an application of the theoretical ideas introduced in this unit (in 
particular, Theorem 1 of Section 3). It concerns a possible approach to the formal definition 
of the trigonometric functions. 


We define the functions S and C to be the solutions of the differential equation 
u"(x) + u(x) =0 (1) 
satisfying the conditions, 
S(0)=0 S’(0)=1 
and 
: c(0)=1 C’(0)=0. 
respectively. 
(i) | Explain why this is sufficient to define the functions S and C. 


(ii) | Prove: (a) that S’ is a solution of Equation (1); (b) that S’(0) = 1 and S”(0) = 0. 
Deduce that S’ = C. 


(iii) Use an argument similar to that in (ii) to prove that C’ = —S. 
(iv) Find the derivative of C? + S?, and deduce that 
C7 +S? =1. 


[Solution on p. 60] 


We could continue to deduce 
other properties of the 
trigonometric functions in the 
same sort of way, but this is 
sufficient to give you the 
flavour of the approach. 
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Appendix 


Solution to the exercise in the Introduction 


1. (i) Equations (a), (c) and (d) are second-order. Equation 
(b) is first-order. 


(ii) Equations (a), (b) and (d) are linear. 
(iii) Equations (b) and (c) have constant coefficients. 
(iv) Equations (a) and (c) are homogeneous. 


Solutions to the exercises in Section 1 


1. The solutions are given in the following table: 


Case 
(Procedure 1.1) 


Auxiliary equation General solution 


y = Ae* + Be™* 

6 = Acos 3t + Bsin 3t 

z = Ae + Be 2 

y = Ae * + Bxe * 

u =e ~(Acos2t + Bsin2t) 
y=A+ Be**? 


> MS MA »AA, 


In each case A and B are arbitrary constants. 


2. (i) The auxiliary equation is 
A* + 2A + w* =0 
ie. (A+) =0. 


This is Case (ii), with the only root 4 = —w. The general 
solution is therefore 


y=e (Ax + B), 


where A and B are arbitrary constants. 


(ii) The auxiliary equation is 
5A? + 60d + Sw? = 0 
thus 


36m? — 100c? 


10 


~6w + 
A= = 


=~(-3 + 4i), 


This is Case (iii) of Procedure 1.1, and the general solution is 
therefore 


3 - . 4 
x = exp (—=at)(Acos —wt + Bsin—at). 
5 5 5 
where A and B are arbitrary constants. 
3. If we can find two linearly independent solutions, then we 


can use Theorem 2 to write down the general solution since, 
with the domain (x > 0), the differential equation is of the 


form specified in the Theorem. To find some solutions, we 


d 
try y = x", as suggested. If y = x” then = = nx"~! and 
2 


= = n(n — 1)x"~?. Thus 


d*y dy 
op Se, ob ie ae -2 __ -1 
2x 72 + + y = 2x* n(n — 1)x" xx” * + & 
= x"(2n(n —1) —n+ 1] 
= (2n* — 3n + 1)x". 
This will equal zero (for all x) if and only if 
es 2n? —3n+1=0. 
roe Ge. iin tha G 


Hence n = $ or 1. So y = x and y = x’? are solutions of the 
given differential equation, and these are linearly 
independent. By Theorem 2 the general solution is 


y=Ax+Bx'/* (x>0) 


4.(i) The auxiliary equation is 3A — 1 — 2A? = 0. 
(ii) This can be rewritten as 

247 — 34 +1=0 
ie: Ol TT 1) oO 


so A=1orl. 


(iii) By Procedure 1.1 the general solution is 
y = Ae + Be** 


where A and B are arbitrary constants. 


5. The solutions are given in the following table: 
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Case 


Auxiliary equation 


A*4+244+2=0 —1+i 
(ii) | A?7-—16=0 +4 
(iii)} A744 +4=0 2 (only) 
(iv)} A427 =0 0 (only) 


In each case A and B are arbitrary constants. 


6. The general solution of the differential equation contains 
sine and cosine terms only in case (iii) of Procedure 1.1, 
when the solutions of the auxiliary equation are complex. 
The auxiliary equation is 


47+ 4kA+4=0 

which has complex solutions when 
16k? — 16 <0 

ie. k? <1 


Le. —1<k< 1. 


7. No. For Theorem 2 to apply we require that the 
d’y 


coefficient of —= never takes the value zero. But x? = 0 when 


dx? 
x = 0, and in this instance x = 0 has not been excluded from 
the domain. 


However we can say that 
(x > 0) 


is the general solution of 


y = Ax + Bx? 


@ d 
2g ~ ee + = 0 (x > 0). 


: dx? d 


Solutions to the exercises in Section 2 


1. (i) The associated homogeneous equation is 

.; + 4y=0. 
Solving this gives the complementary function 

y = Acos 2x + Bsin 2x. 
(ii) Trying a solution of the form y= C_ (C constant) in 
=i + 4y = 8 gives 0 + 4C =8, that is C = 2. Thus a 
particular solution is 

y =2. 


(iti) By Theorem 1 the general solution (of the 
inhomogeneous equation) is 


y = Acos 2x + Bsin 2x+2. 


2. (i) For y = mx + n to be a particular solution, we require 
(for all x). 


This can only be true for all values of x if m and n satisfy the 
simultaneous equations 


—2m + 2(mx + n) = 2x + 3 


—2m + 2n = 3 (equating constant terms) 
2m = 2 (equating coefficients of x). 
Solving for m and n we obtain 
5 
m=t1,n 


5 


(Procedure 1.1) 


General solution 


y=e *(Acosx + Bsinx) 
y = Ae* + Be ™ 

y = Ae™ + Bxe?* 
0=A+Bt 


So a particular solution is 
ered 
zn 
(ii) For y = mx + n to be a particular solution, we require 
2m + (mx + n) = 2x. 


This can only be true for all values of x if m and n satisfy the 
simultaneous equations 


2m+n=0 (equating constant terms) 

m = 2 (equating coefficients of x). 
Substituting m = 2 into 2m + n = 0 gives n = —4,soa 
particular solution is 

y=2x -4. 
3. If we try 
y= me", 
d d? 
then oe —me ~*~ and ss = me ~*~. So we need 
dx dx 


20° == 2me-* + 2me” * + me. 


2 
This is satisfied if 2 = 5m, that is m = 5 So a particular 


solution is 


=-e *. 
3 


4. We can either try a sinusoidal solution directly or use the 
phasor method. 


Method 1: We try 


y =mcos 3t + nsin 3t. 


then 
d 
ae. = 3ncos 3t — 3msin 3t, 
dt 
and 
d2 
=, = —9mcos 3t — 9nsin 3t. 


SO we require 
(—9m + 9n)cos 3t + (—9n — 9m)sin 3t = 9 cos 3t. 


This can only be true for all values of t if m and n satisfy the 
simultaneous equations 


—9m + 9n =9 
—9m — 9n = 0. 
Solving these equations gives m = —4, n = }. So the 


particular solution is 


y = —$c0s 3t + 4sin 32. 
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Method 2: We rewrite the right-hand side of the differential 
equation as Re(9e*") and try y = Re(ze*"). Using Procedure 
2.2(a) we obtain 


(—9 + 9i)z = 9, 


) 1 
so Z = ——__ = ——__ = — 
—9 + 9i 


is therefore 


eee 5 - 5 The particular solution 


y = —4cos 3t + 4sin 3. 


5. Method 1: We try 
y = mcos2t + nsin 2t. 
Then 


d 
© Ines — Pan 2 


dt 
and 
dt? 


So we need 


= —4mcos2t — 4nsin 2t. 


— 8mcos 2t — 8nsin 2t = 8 cos 2t + 16sin 2t. 


This can only be true for all values of t if m = —1 and 
n= —2. Thus a particular solution is 


y = —cos 2t — 2s1n 2t. 


Method 2: We write the right-hand side of the differential 
equation as Re((8 — 16i)e7") and try y = Re(ze?"*). Using 
Procedure 2.2(a) we obtain 

(—4 — 4)z = (8 — 16i). 
So z = —1 + 2i and the particular solution is 

y = —cos 2t — 2sin 2t. 


6.(i) We can see that y = cos 2t is a solution of 
2 


dy 
dt re % 
alternative methods of solution are given below. 


Method 1: We try 


+ 4y = 0, so we have an exceptional case. Two 


y = t(mcos 2t + nsin 2t) 
then 


2 = (mcos 2t + nsin 2t) + t(2ncos 2t — 2msin2t) 
and 
2 


= = (4ncos 2t — 4msin 2t) + t(—4mcos 2t — 4nsin 21). 


2 
Substituting into = + 4y gives 4ncos 2t — 4m sin 21, so we 
require 
4ncos 2t — 4msin 2t = 2 cos 2t. 


This can only be satisfied for all values of t if m = 0 and 
= 4. So the particular solution is 


y = 4tsin 2t. 


Method 2: We rewrite the right-hand side of the equation as 
Re(2e*") and try 


y = Re(tze”"*) 
then 
dy 


a = Re((l + 2it)ze*") 
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and 
2 


d 
— — Re((—4t + 4i)ze2it) 
dt 

2 


me ied 
Substituting into —> + 4y gives Re(4ize7"), so we need 


dt ek: 
4iz = 2. 
Thus z = —4i and the particular solution is 
y = 3tsin 2t. 


(ii) 4 = 1 is a solution of the auxiliary equation 

4* — 34 +2 =0. So y = e* is a solution of the associated 
homogeneous equation, and we have an exceptional case. 
Following Procedure 2.3 we try 


y=mxe". 
| d d? 
Then a m(1 + x)e* and oes = m(2 + x)e* and so 
dx dx 
d’y 


dy 
23-3. t a= (2 + x) 
Thus we have a particular solution if m = 
solution is therefore 


y= —4e. 


— 3(1 + x) + 2x)me* = —me*. 


— 4. The particular 


7. In this case the general solution of the associated 
homogeneous equation is y = Ae* + Bxe*. So not only is 

y = me* going to reduce to zero if substituted, but so is 

y = mxe“. So we need to try something else to find a 
particular solution. The trick is to multiply by x again. That 
is, we try 


y = mx*e’. 
dy a d*y 
Then te = (2x + x*)me* and Za = (2x + x? +2 + 2x)me*. 
So 
dy dy ; 
ear See, Le oe 4 _ 1,2 2 
7 at (x? + 4x + 2 — 2(x? + 2x) + x?)me* 


= 2me. 


Thus we have a particular solution if m = 4. The particular 
solution is therefore 


y = ix*e. 
8.(i) We know from Theorem 2 that we can deal with the 
terms on the right-hand side of the differential equation 
separately, and we know from Exercise 7 that y = 4x?e* is a 


particular solution for the case where the right-hand side is 
e*. 


Let us now deal with the e?* term. To find a particular 
solution of 


d 
a7 ~ 2a ty = e%, (1) 
x 


we can try y = me**. This satisfies the equation if 
(4m — 4m + m)e?* = e? 
That is, if m = 1. So a particular solution of Equation (1) is 
y= 
Using Theorem 2 a particular solution of the given 
differential equation is, therefore, 


= 4(ix7e*) — 3e?* = 2x7e* — 3e?*. 


52 


(ii) We can find a solution by trying 
y=1+mcos4t + nsin4t 
(or by the complex number method). 
This is a solution if (substituting in the differential equation) 


4(—16mcos 4t — 16nsin 4t) + 9(—4msin 4t + 4ncos 4t) 
+ 100(/ + mcos 4t + nsin 4t) = 9cos4t + 100. 


That is, if 
100] + (36m + 36n) cos 4t + (— 36m + 36n) sin 4t 


= 9cos 4t + 100. 
This can only be true for all values of t if _ 
36m + 36n = 9 (equating coefficients of cos) 
— 36m + 36n = 0 (equating coefficients of sin) 
100/ = 100 (equating constant terms). 


Solving these simultaneous equations gives | = 1, 
n= and m=#. Soa particular solution is 


y=1+(cos4t + sin 4t). 


9. (i) First find the complementary function, by solving 
a6 
—~ + 40 =0. 
we 
This has auxiliary equation 
77+4=0 


which has solutions 4 = +2i. Hence the complementary 
function is 


6 = Acos2t + Bsin2t 
where A and B are arbitrary constants. 
To find a particular solution, try 
6 =mt +n. 
This satisfies the given inhomogeneous differential equation if 
(for all t); 


that is, ifm = 4 and n=0. 


4(mt + n) = 2t 


So the general solution of the given equation is 


6 = Acos2t + Bsin2t + 4t. 


(ii) The associated homogeneous equation is 


F2 
a7 -¥=0 
Its auxiliary equation is 
347 —2A —1=0 
which has solutions 4 = 1 and 4 = —}4. Hence the 


complementary function is 
Y= Ae* + Be ®, 
where A and B are arbitrary constants. 
To find a particular solution, try 
Y= me** +n. 
This satisfies the given equation if 
3(4me?*) — 2(2me?*) — (me?* + n) = e?* + 3. 
That is, if m =4 and n= —3. 
So the general solution of the given equation is 


Y= Ae i, + 5e?* — 3. 
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(iii) The auxiliary equation of the associated homogeneous 
equation is 


A*-344+2=0 


which has solutions 4 = 1 and 4 = 2. So the complementary 
function is 


y = Ae* + Be”, 


where A and B are arbitrary constants. Note that in finding 
a particular solution we have an exceptional case. Both of 
the functions on the right of the given differential equation 
are solutions of the associated homogeneous equation, so for 
a particular solution we should try 


y =mxe* + nxe”*, 
Then 


d 
= = (1 + x)me* + (1 + 2x)ne”* 


and 
2y 
re = (2+ x)me* + (4+ 4x)ne?*. 

x 


. . . dy dy - 2x 
Substituting into —; — 3— + 2y gives —me* + ne™ so we 
dx dx 
require 
—me* + ne* = 2e* — Se". 


This can only be true for all values of x ifm = —2 and 
n = —5. Thus the general solution of the given differential 
equation is 


y = Ae + Be** — 2xe* — 5xe*. 
(iv) The auxiliary equation of the associated homogeneous 
equation is 

A? + 2wd + 2m? = 0. 
Thus 

= 1(—2M + \/4w? — 8w”) = w(—1 + i). 

So the complementary function is 

y=e “(Acosat + Bsinot). 


I will deal with the two parts of the right-hand side of the 
differential equation separately. (This is not necessary—I just 
prefer this way here.) A solution of 

d*y dy 
72 - 207 + 2w*y = wh 


can be found by putting y = m (where m is a constant), 
which yields the particular solution 


y = gh. 


A solution of 


d? d 
ay + Ww + 2w*y = Ssin 2a, (1) 


can be found by putting 
y =mcos 2wt + nsin 2mt 


(or by complex numbers). Then 


d 
~ = 2wncos 2wt.— 2mmsin 2mt 


—— = —4qw*mcos 2ut — 4w?nsin 2at. 
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Substituting into the left-hand side of the differential 
equation gives 


(—2m + 4n)w? cos 2wt + (—4m — 2n)w? sin 2ot. 


We require this to equal 5sin 2wt, which means that m and n 
must satisfy the simultaneous equations 


—2m+ 4n=0 
—4m — 2n = 5/w’. 
; , 1 1 
Solving these equations gives m= —— andn= — ==, BO 
7) 2m 


a particular solution of Equation (1) is 
1 
y = — =, (2cos 2wt + sin 2wt). 
2w 


Using Theorem 2, a particular solution of the given 
differential equation is 


1 
y= — Fpp2 (2 cos 2ot + sin 2ot) + th. 


Hence the general solution is 


1 
y=e “(Acoswt + Bsinwt) — 552 2 cos 2mt + sin2at) + th. 


10 (i) y = me*. 
(ii) y = mcos 3x + nsin3x (or y = Re(ze***)). 


(iii) Note that the complementary function is 
y = Ae** + Be ?*, so we need to try y = mxe~?* for a 
particular solution. 


(iv) Note that the complementary function is 
y = Acos2x + Bsin2x, so we must try 


y =kxsin2x + Ixcos2x +mx+n 


for a particular solution (or y = Re(zxe?'*) + mx + n). 


11 (i) The complementary function is 
y=e *(Acosx + Bsinx). 
A particular solution is 
y=2., 
So the general solution is 
y=e *(Acosx + Bsinx) + 2. 
(ii) The complementary function is 
y = Ae** + Be *. 


For a particular solution, we try y = Re(ze**). By Procedure 
2.2(a) 
(-—9-—4)z= -i 
Le. Z=—i. 
13 
Thus a particular solution is 


1 
= ——sin 3x, 
a 
Hence the general solution of the given differential equation 


1S 


1 
y= Ae2* + Be~2* =A 0m 3x. 


Solutions to the exercises in Section 3 


1. (i) Do not be thrown by the different notation I have used 
here. The question could equally well be written; 

d’y dy Tt 
ga to p= Cee = ate 
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where y = u(t). (Letters other than y and t can be used for 
the variables, of course.) The auxiliary equation is 


47+4=0. 
Therefore 4 = +2i and the general solution is 
u(t) = Acos2t + Bsin 2t. 
Now u(z/2) = 0 and so 
0 = Asinn + Bcosn 
=Ax0+Bx -l. 
Therefore B = 0 and hence 
u(t) = Asin 2t. 
Now u'(t) = 2A cos 2t and so the condition u’(x/2) = 1 gives 
1 = 2Acosa = —2A. 
ie. A = —3. 
The required particular solution is therefore 


u(t) = —4sin 2¢. 


(ii) The complementary function is 
y = Ae + Be** 
(see Frame 1 for details). 
A particular solution is 
ye? 
The general solution is, therefore 
y = A& + Be** + 2. 
When x = 0, y=4 and so4=A+B + 2. That is, 


A+B=2. (1) 
d 
Now dy = Ae* + 2Be** and when x = 0, Peis —1. Thus 
dx dx 
—1=A+2B. (2) 
Solving the Simultaneous Equations (1) and (2) gives A = 5, 
B = —3. The required particular solution is therefore 


y = 5e* — 3e* + 2. 
2. This is covered in Frame 7 of the audio-tape. 


3. Problems (i) and (iii) are initial condition problems, (ii) 
and (iv) are boundary condition problems. 


In each case, the differential equation is the same. Its general 
solution is 


u(x) = Acos2x + Bsin 2x. 
The derivative is 
u'(x) = —2Asin2x + 2Bcos 2x. 
Ni solutions satisfying the given conditions are found 
Ow. 


(i) The condition u(0) = 1 gives A = 1. The condition 
u'(0) = 0 gives 0 = 2B. That is B = 0. The required solution 
is therefore 


u(x) = cos 2x, 


(ii) The condition u(0) = 0 gives 


0=A+0. 
The condition u(z/2) = 0 gives 
0=A-Q0O. 
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Thus B is arbitrary and A = 0. Any solution of the form 
u(x) = Bsin 2x 
satisfies the problem. 


(iii) The condition u(0) = 0 gives A = 0. The condition 
u'(0) = 0 gives B = 0. So the required solution is 


u(x) = 0. 


(iv) The condition u(—2) = 1 gives A = 1. The condition 
u(x/4) = 2 gives B = 2. So the required solution is 


u(x) = cos 2x + 2sin 2x. 


4. (i) The required solution is 
y= 
as we have a homogeneous differential equation with initial 


d 
conditions of the form y = 0 and = == (}. 


(ii) The complementary function is 
y= Ae + Be ™*. 
A particular solution is 
y= -8. 
So the general solution is 
y = Ae + Be * — 8. 


The condition.‘y = 0 at x = 0’ gives 


A+B-—8=0. - (1) 
dy = ae... 
Now — = Ae* — Be *, and so the condition — = 0 at 
dx dx 
x = 0’ gives 
A—B=0. (2) 


Solving the Simultaneous Equations (1) and (2) gives 
A = B= 4. Therefore, the required particular solution is 


y = 4e* + 4e°* — 8. 
(iii) The general solution is 
x = (At + Bye ™ 


(see Exercise 2(i) of Section 1). The condition ‘x = aatt=0 
gives 


a=B. 


dx 
Now Pier w(At + Bje'™ + Ae ™, and so the condition 


e 


a 0 att =O gives 


0= —wB +A. 


Thus A = wB = aw. The required particular solution is 
therefore 


x = a(wt + lje™™. 


5. (i) For first-order equations the theorem is as follows. Let 


d 
p(x) + a(x)y = fle) 


be any linear first-order equation for which p(x) does not 
take the value zero for any value of x. Then there is one and 
only one solution of this differential equation which also 
satisfies a condition of the form y = b when x = a. 
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We can find some soluutions to the differential equation in 
(ii) by the following incautious manipulation. Separation of 
variables gives: 


[? = (e 
y Ix 
so log.y =log.x+C 


ie. p= Ax: 


I say incautious because the manipulation performed is 
invalid if x or y is zero or negative. However it happens that 
the solutions found here are alright for all values of x and y, 
as I can show by checking. 


Check: If 
y = Ax 
then 


for all real x, as required. Furthermore, if x = 0, then 
y = 0 (whatever A is). Hence all functions of the form 


y = Ax 


satisfy the given differential equation and initial condition. 


(iii) There is not a unique solution in this case. We do have a 
linear equation here. What is wrong is that the condition 
concerning p(x) is not satisfied (for p(x) = x is zero when 

x = 0). This example illustrates the fact that the condition 
p(x) # 0 in the theorem is crucial. Without it, the theorem 
does not work. 


6.(i) The general solution is 
u(t) = e (A cost + Bsint). 

The condition u(0) = 0 gives 0 = A and hence 
u(t) = Be *sint. 

Thus- 
u'(t) = Be~~ (cost — 2sint), 


and so the condition u’(0) = 2 gives 2 = B. The required 
solution is therefore 


u(t) = 2e~ *' sin t. 


This is an initial condition problem. 


(ii) The general solution is 
u(t) = Acos 3t + Bsin 3t. 

The condition u(0) = 0 gives 0 = A and hence 
u(t) = Bsin 3t. 

Thus 
u’(t) = 3Bcos 3t 


and so the condition u’(z/3) = 1 gives 1 = 3Bcosz = —3B. 
That is B = —4. The required solution is therefore 


1 
u(t) = —, sin 3t. 


This is a boundary condition problem. 
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(iii) We have a homogeneous equation, with initial 
conditions equal to zero. The required solution is therefore 


y = 0. 


This is an initial condition problem. 


Solutions to the exercises in Section 4 


1. The phasor of cos 5t is 1 and that of sin 5t is —i, thus the 
required phasor is 


z=4+4 3i. 
The amplitude is 


lz| = /47 + 3? =5. 


The phase is 


Arg z = arctan (2] = 0.64. 


2. The auxiliary equation is 
A? + 2a, + wm? =0 
which has roots 


A= —a+t./o? — w?. 


The roots will be real provided a? — w? > 0. 


3. The graph is a sinusoid curve of angular frequency w (the 
amplitude and phase depend on the initial conditions). 


4. (i) We require «* — w” to be negative and « to be positive 
(see table on page 36). Thus we require 


O<a<a. 


(ii) (a) For this equation « = 1.5 and w = 3. Thus in this case 
0 < a <q, and so the solution is a decaying oscillation. To 
ensure that the solution also fits the given initial conditions it 
must take the form shown in Figure 1 (starting at y = 4 with 
zero gradient). 


~ Initial gradient= 0 


Figure 1 


(b) For this equation « = 4 and w = 2. Thus in this case 

a > w > 0 and so the solution is a combination of decreasing 
exponentials. To satisfy the given initial conditions the 
solution must take the form shown in Figure 2. 


Initial gradient = | 


Figure 2 
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5. The solution is the sum of (J) a particular solution, which 
will be a sinusoid, and (I/) the general solution of 


d? d 
aa + 2a + wy =0. 


Since a > 0, this general solution will be decreasing, and will. 
tend to zero as x becomes large and so is a transient. Since 
a <w this transient is oscillatory. The initial conditions only 
affect the transient term in the solution. So whatever the 
initial conditions the long-term behaviour of the solution is 
the same. It is a sinusoidal oscillation, of frequency v. The 
amplitude and phase of this oscillation depend on v, « and w. 
(Example 1 on page 39 illustrates how the amplitude and 
phase can be calculated.) 
6. The general solution of the differential equation is 

y = Acos 3x + Bsin 3x. 


The angular frequency of this oscillation is 3. For the 
condition ‘y = 4 when x = 0’ to be satisfied we require 


d 
A= 4. Now = = —3Asin3x + 3Bcos3x and so for the 


condition 2 = 9 when x = 0’ to be satisfied we require 3B 
= 9, that is B = 3. Thus the required particular solution is 
y = 4cos 3x + 3sin 3x. 
The phasor (see Unit 5) of this oscillation is 
z=4—- 3i. 
Hence the amplitude of the oscillation is 
lz] = ./47 + 37 =5, 
and the phase is 
Argz = —arctan (3) = —0.64. 
7. All the solutions of u”(x) + w?u(x) = 0 are sinusoidal 
oscillations of angular frequency w. These solutions repeat 
themselves after a period = and so any solution u(x) of this 


differential equation satisfies the condition 
2n 
u(x) = uj x + — 
@ 
for all values of x. In particular, when x = a we have 


| 2n 
u(a) = ula + — 
7) 


A solution can therefore never be found satisfying boundary 
conditions of the form 


: 2 
u(a) =k and ufo + = == | 
w 
if k #1. On the other hand, if k = |, then any solutions 
2 
satisfying u(a) = k will also satisfy u f + =| =1(=k). 


Consequently the two conditions are equivalent in this case, 
and so will only place one condition on the two arbitrary 
constants in the general solution. So for k = I, the solution is 
not unique. 


8. From Example 7 of Section 1, the general solution is 


x = Ael-2+V3)a + Bel 2-V3)ar 


For this to satisfy the initial condition ‘x = a at t= 0’ we 
require 


a=A+B. (1) 


Now 


d 
= = (-2+ /3)mAe- 2+ — (2 + ./3)mBe 2-9 


6 


For this to satisfy the initial condition eile 0 att = 0 we 


require 


0 = (./3 — 2)wA — (2+ /3)mB. (2) 


Solving (1) and (2) gives 


(./3 —2)A = (2 + ./3)B 
= (2 +/3)(a — A). 
Thus 


a oye -(} Bla 


ofthe fh 
So the required solution is 


roe oe. 
— alin MIN An 2t Vd 4 [EV A A - 203 ee 
a [+ 3 r + F 3 r } 


and 


9. The solutions of the auxiliary equation 
A? + 2BwdA + w*? = 0 


are 


—2Ba + ./4B?w? — 40? 


Z 
= w(—B + \/p? — 1). 


Hence (since f? < 1) the general solution is 


x = e814 cos (w./1 — B?)t + Bsin(w./1 — B?)t 
For this to satisfy the conditions ‘x = a at t = 0’ we require 
A=4a4. 
Now 


d 
vane — Bwe 


dt 
(Acos (w./1 — B?)t + Bsin(w./1 — B?)t) 
(—Asin(w./1 — B7)t + Bcos(w./1 — B?)t). 


— pot x 


6 


ae ; 
For this to satisfy the condition a= 0 at t = 0 we require 


w(—AB + B/1 — B?) =0 
ie: 
2 a7 Ap ap 
es Se es 2 


So the required particular solution is 


x = ae~** {cos (w./1 — B?)t + ————=sin (w,/1 — B?)t}. 


TF 


This is a decaying oscillation, as sketched in Figure 3. 
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Initial gradient= 0 


Figure 3 
10. To find the phasor z of the steady-state solution, we try 
y = Re(ze>"). This gives 

(4(5i)? + 9(5i) + 100)z = 9 


9(Si)z = 9. 
So 
1 es 


z=—=-- 


Thus the amplitude of the steady-state solution is 
1 
F4 a 


5 


and the phase is 
A us 
rgz= —~. 
ro ed 

Its angular frequency is 5. [The actual steady-state solution 


1 Tl 
is, therefore, x = — St ——}. 
is, therefore, x 500s | 


11. If we try y = Re(ze™) we obtain the following equation 
for the phasor z of the steady-state solution: 


(—v? + 2Bo(iv) + w)z = 1 


2 1 
~ (w? — v7) + (2Bav)i- 
The amplitude of the steady-state solution is 
1 
ica Ph a. eee, FRY 
Jo — v*)” + (2Bov) 


12. The nature of the solutions of the differential equation 
d’y dy 
+ 4 +@ 0 
re Cag a 
depends on the nature of the roots of its auxiliary equation, 
A? + 2ad + w? = 0. 
These are 
A= et fo =o’. 
They are real if «? — w? is positive, complex if «7 — w? is 
negative. Let us look at the various cases in turn. 
(I) « positive, «? — w? positive 


The roots are real, and both are negative. (The root 
—a + ./a? — w? must be negative because « > ./a? — w). 
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Thus the solution of the differential equation is a linear 3. (i) 
combination of e*'* and e*”* with 4, and 4, negative. Thus 
the typical solution in this case is a decreasing exponential. 


(II) « negative, a? — w? positive 

Reasoning similar to that in Case (I ) shows that in this case 
the solution is a linear combination of e*'* and e*”* with A, 

and A, positive. Thus the typical solution of the differential 

equation in this case is an increasing exponential. 


(III) « positive, «* — w? negative 
In this case the roots of the auxiliary equation are complex: 


A= —ati,/w? — a’. 


The solutions of the differential equation are linear 
combinations of e~** cos (./@? — a? x) and 

e~ sin (./@* — a x). The typical solution is decreasing and 
oscillatory. 


(IV) « negative, «? — w* negative 

The formulae for the solutions are as in Case (III), but now 
e * is an increasing exponential. Thus the solutions of the 
differential equation are oscillations of increasing amplitude. 


The Euler approximation to y(3.2) is shown in bold type. 


2 


d 
2 ey ia 0 is 


(ii) The general solution of de 


; f ; : y = Acosx + Bsinx. 
Solutions to the exercises in Section 5 poe Lo 5) 
dy The condition ‘y = 0 when x = 0’ gives A = 0. The condition 
1. (i) — =z ‘dy 
dx 
dx 


= 1 when x = 0 gives B = 1. So the required particular 


ae = —4z — Sy + cosx solution is: 
dx 

- Sy y =sinx. 

ti) . (iii) The results are shown graphically in Figure 4. The 
te general shape of the results is correct—a half cycle of an 
ee es 32 sin y. oscillation. However the results are not precisely accurate. 


(The values of y should not rise above 1). This inaccuracy is 
hardly surprising with so long a step-length as h = 0.2. 


d 
2. In each case, replace = by (y,41 —y,)/h, 


dz 
in by (z,-1; — Z,)/h, and x, y and z by x,, y, and z,. HH PL iT iat HEE zs 


SEBSRASRSY CRSA KR ESE SRA Se 
(i) For these equations we obtain “By H 


yr+1 — Yr eo 
ee eee, 
and 
on+i Zz, 
h = —4z, — 5y, + Cosx,, SoBe SERS ee a\a 
é s seGeeneunecuaan 
SPRRE ESRRPRERE RSME eee 


ESSER ERERELERERL 
ao I 


which can be rearranged to give the pair of recurrence SSSSSnSnanERN 


relations 


Yr+1 ae: yp az hz, 
2,41 = —dShy, + (1 — 4h)z, + hoosx,. 


(ii) In this case we obtain 


Vr+1 ee i 


—2Z 
<= —32siny,, 


which can be rearranged to give the recurrence relations 


Vr+1 ay Yr hz, 
2,41 = 2, — 32hsiny,. 
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The Euler approximation to y(1) is shown in bold type. The 
true solution is given by 


1 
y= g(lles™ — 3e7 10), 


Figure 5 shows a comparison of this true solution with the 
Euler approximations. 


SESCSESRESTNRHTRESRe sande CeeTeSaaks Bases 
ECP EEE EES 
de stazeet tatasee eta easeeeete 


oe 
EO 


> 


HE 
Ht 
seees 


Figure 5 


The method has failed miserably to give results that are 
anywhere near reasonable. However, we need not despair 
since if we use a small enough step-length, we can get useful 
results. Clearly, we have to be careful when we use numerical 
methods. 


1.6z, — 0.4y,+ 0.4 ° 


The Euler approximation to y(1) is shown in bold type. The 
true solution is given by 


y=e™—e +4, 


The graph indicates that the numerical results have the same 
shape as the solution curve. 


a ae eetaaeeeeia 


i: sc nail 


ie Hee ta 
i een EEE EHEEHE: 


gle as Sueeenuene 


p Se 
2818 RECRSSRewS S 
sf jessesssiy’ 


; ae oe Soe Oo = 


wrewipaces 


Figure 6 
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Solutions to the exercises in Section 6 

1. (i) The auxiliary equation 
47 +614+5=0 

has roots 4 = —5, 4 = —1 and so the general solution is 
yes Aer 4 Ber, 


The initial condition ‘y = 2 at x = 0’ is satisfied if 


2 =A + B. (1) 
ee eee 
The condition oe 1 at x = 0 is satisfied if 
1 = —5A —B. (2) 
3 11 
Solving Equations (1) and (2) gives A = 2 a rt The 


required particular solution is therefore 
1 
= —(lle-* — 3e7 >*); 
yaa ) 


(ii) The complementary function is 
yoo = + F 
A particular solution is y = me?*, where 


4me** + 4me?* = e?* 


1 
ie. Bi= 3 The general solution is therefore 
~- l 2x 
y=Ae *+B+ g° ° 


d 
The initial conditions ‘y = 0 and = =latx=0O give 
x 


1 
0O=A+B+ 8 
1= —-2A+ 2 
— ; 
3 
Solving these simultaneous equations we obtain A = a 


1 
B= rr The required particular solution is therefore 


a | ee 
oe a 


(ili) The complementary function is 
y = Ae* + Be7 *. 


A particular solution of 


d’y 
— — 9y = 18 
dt’ : 
is 
y= —2, 


A particular solution of 
d2 
ie ~ 9y = 3sin 31 

is a sinusoid with phasor z, where (putting y = Re(ze>")) 
—18z = —3i, 


; 3 
Le. z = —i. 
6 


That is 
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1 
= —-—sin 3t. 
an 
The general solution is therefore 


y = Ae* + Bo* —2— g sin 3 


d 
The conditions ‘y = 0 and & = 0 att=0 give 
0O=A+B-2 
1 
; : 13 11 
Solving these simultaneous equations gives A = rr B= re 


The required particular solution is therefore 
| 
y= 7pl3e" + 1le~** — 24 — 2sin 3¢). 
(iv) This can be solved by direct integration. We have 


d 
aaa ela iag 
u 


d 
We requir - = 2 when u = 1, so A = 3. That is 


dé 1 
Se 3 -— 
Hence 
d= ff us au 
u 
= 3u — log .u + B. 


The condition ‘6 = 1 when u = 1’ gives 
1=3+B. 

Thus B = —2 and the required particular solution is 
0 = 3u — log .u — 2. 

2. (i) The complementary function is 
y = Acos4t + Bsin 4t. 


This is an exceptional case. We find a particular solution (of 


the form y = t x (a sinusoid)) by putting y = Re(zte**). Then 


yn ; 
at = Zé (1 + 4it), 
and 


2 


d 
at = ze*#(4i + 4i(1 + 4it)) = ze**(8i — 160). 


So 
a 
dt? 

Thus the phasor z satisfies 


8iz = 4 — Bi, 


+ 16y = 8ize*#. 


1. ; Sg 
Le. z= —1 — 5 So a particular solution is 


1 
y = t(—cos4t + 5 sin 41). 


Hence the general solution is 


1 
y = Acos4t + Bsin4t — tcos4t + at sin 4t. 
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(ii) The complementary function is 
y= Ae “ + Bie”. 

To find a particular solution, try 
y=It?e"* + mt+n. 

Then 


dy 2),- 28 
am et — 2t je“ + m 


and 


PA 
ed = (2 — 4t — 4t + 4t?)e~7* = (2 — 8t + 4t2)e~ 2". 


Thus 


d d | 
ar mn 4 + 4y =le~**((2 — 8¢ + 41?) + 4(2¢ — 22) + 42 
+ 4m + 4(mt + n) 
= 2le~** + 4mt + 4(m + n). 
So we require 


2le~** + 4mt + 4(m + n) = 3e77* + 8t. 


Therefore 2/ = 3, 4m = 8, 4(m + n) = 0. That is 
[= - m = 2, n= —2. 


The general solution is therefore 


3 
y = Ae~* + Bte~** + ae + 2t — 2. 


3. (i) Since (67 -4x1x1)>0 and 6>0 the graph 
has the form 


y 


Initial gradient =—2 


Figure 7 


(ii) Since (17 — 4 x 6 x 1) <0 and 1 > 0 the graph has the 
form 


Initial gradient = | 


Figure 8 


d 
4. (i) Since the coefficient of _ is positive the 


complementary function is a transient. The long-term motion 
is therefore independent of the initial conditions. The non- 
transient term will be a particular solution of the form 


x = Acos(3t+ @)+m 


where A, ¢ and m are constants. It is easiest to deal with the 
constant term m separately. For x = m to be a solution of 
d*x dx 
ge oe ae 2 
dt? edt 


1 
we require 4m = 2. That is m = 5 and the solution is 


ae 1 
<. 
To. find the sinusoidal solution of 
d*x dx 
Ph oe 5, + 4x = 2sin 3t, 


we can use phasors. The phasor z of this solution is found by 
putting x = Re(ze*") (remembering that the phasor of 2 sin 3t 
is —2i). So 


(3)? + 3i + 4)z = —2i 
—2i 25 +31) 2 


16. 2 = ors 34 343 + 58). 
Thus A (the amplitude) is given by 
A = |z| = 2 


Ta 


and ¢@ (the phase) is given by 
5 
o@ = Argz=n — arctan, = 2.11. 


Thus the steady-state solution is 


x= ; + eee (3t + 2.11). 


34 


The graph is shown in Figure 9. 


Figure 9 


The smallest value taken by x during this part of the motion 
is 


1 z 


Oe. 


This is positive, so x does not take negative values. 


(ii) ‘a long time after the object is set in motion’ means 
‘sufficiently long after the motion has started for the transient 
term to have become negligible’. We can spell this out more 
precisely if we find the transient. This is the complementary 


function, which is 
“fie ae 15, 


— ptt 
x=e A cos —t + Bsin ~— 
2 2 


The term that determines how rapidly this transient dies 
away is clearly, e~ *'. Precisely when it is reasonable to 
neglect this term depends on the particular problem under 
consideration. However, it would often be adequate to 
require that e~ ** < 10~*. This occurs when t > 2 log, 10+ 
= 18.4. 
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5. The general solution of the differential equation is 
u(x) = AcosAnx + Bsin Anx. 


For u(0) = 0 we need 


0=A. 
For u(1) = 0 we then need 
0 = Bsin An. 


For a non-zero solution we require that B is not zero, so we 
need 


sin Ax = 0. 


This occurs if A is an integer. 


6. (i) This is a consequence of Theorem 1 of Section 3. The 
differential equation u”(x) + u(x) = 0 is of the form specified 
in that theorem, so there is exactly one solution of this 
equation satisfying a given set of initial conditions. Hence 
there is one and only one function S satisfying the given 
definition, and so this is an adequate definition. A similar 
comment applies to the definition of C. 


(ii) (a) Since S is a solution of Equation (1): 


S”"+S=0. 
Differentiating this gives 
— + Ss’ _ 0. 


Le. (S’)” + S'=0. 


So S’ is a solution of Equation (1). 


(b) Also 

S’'(0) =1 (by the definition of S) 
and 

S”(0) = —S(O) (from Equation (1)) 


== () (by the definition of S). 


So the function S’ is a solution of Equation (1), that satisfies 
the initial conditions 


s‘0)=1, = (S’)'(0) =0. 
But the function C satisfies the same initial conditions, and 


by Theorem 1 of Section 3 there is only one solution of 
Equation (1) satisfying these initial conditions. Hence 


C= S’. 


(iii) C’ is a solution of Equation (1), by an argument just 
like that in part (ii). Also 


CO) =¢ (by the definition of C) 
and 
C”(0) = —C(0) (from Equation (1)) 
= —1 (by the definition of C) 


Now the function —S is also a solution of Equation (1), and 
satisfies the initial conditions 


— S(O) = 0, —$’(0) = -1. 
Hence, by Theorem 1 of Section 3: 


(iv) The derivative of C? + S? is 
2CC’ + 28S’ = 2C(—S) +2SC =0. 
Hence C? + S” must be a constant function. So 


C? + S? = C*(0) + S?(0) 
= 1. 
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